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ABSTRACT 

The  Quantum  1/f  Theory  has  been  developed  under  the  present  Grant 
and  has  been  applied  to  the  noise  suppression  in  various  electronic 
devices.  The  Quantum  1/f  Theory  derives  fundamental  quantum 
fluctuations  present  in  the  elementary  processes  of  physics  at  the  level 
of  the  quantum  mechanical  cross  sections  and  process  rates.  New 
developments  described  in  this  report  include  the  derivation  of  the 
conventional  Quantum  1/f  Effect  in  second  quantization,  and  for  an 
arbitrary  number  of  particles  in  the  final  state.  They  also  include  the 
derivation  of  the  quantum  1/f  cross-correlations  between  scattering 
cross  section  fluctuations  at  different  angles.  These  cross-correlations 
are  needed  for  a  more  exact  calculation  of  quantum  1/f  noise  in  various 
kinetic  coefficients,  such  as  the  mobility  of  the  current  carriers.  Another 
new  development  is  the  derivation  of  a  formula  which  describes  the  effect 
of  the  finite  mean  free  path  of  the  current  carriers  in  condensed  matter  on 
the  quantum  1/f  noise.  Finally,  an  interpolation  formula  linking 
conventional  and  coherent  quantum  1/f  noise  was  developed.  Practical 
applications  to  semiconductor  samples,  pn  junctions,  junction  and  MIS 
detectors,  Josephson  junctions  and  SQUIDs  were  developed  and  are 
presented.  Optimal  design  principles  based  on  the  Quantum  1/f  Theory 
have  been  developed  and  are  described,  applied  and  explained. 
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I.  INTRODUCTION 

Under  the  present  Grant  major  progress  was  achieved  by  the  author 
both  in  developing  the  Quantum  1/f  Theory1'13,  and  in  applying  it  to  high- 
technology  systems16*20,  in  particular  to  low-noise  design  of  electronic 
devices,  including  ultrasmall  devices.  The  present  report  reviews  this 
progress,  presents  the  physical  meaning  of  the  Quantum  1/f  Effect,  and 
discusses  the  impact  it  has  both  on  practical  hi-tech  applications  and  on 
the  foundations  of  quantum  physics. 

The  Quantum  1/f  Effect  was  first  derived  formally  by  the  author  in 
second  quantization  for  a  pair  of  bosons  or  fermions  in  1986.  In  1987  this 
second  quantization  treatment  was  extended  by  the  author  to  an  arbitrary 
number  N  of  bosons  or  fermions.  This  allowed  for  the  first  time  a 
rigorous  derivation  of  the  1/N  factor  present  in  the  Quantum  1/f  Theory. 
Before  1986,  the  simplified  physical  derivations1*2  of  1975  and  1980 
were  used,  similar  to  the  elementary  theory  of  diffraction,  translated  into 
the  time  domain.  The  present  second-quantized  formulation  of  the  theory 
is  physically  equivalent  to  the  previous  formulation,  but  is  more  rigorous, 
while  still  easily  accessible.  The  present  formulation  justifies  the  use  of 
four  single-particle  wave  functions  in  the  autocorrelation  functions, 
provides  a  theoretical  basis  f  me  simplified  elementary  model  based  on 
interference  beats  between  the  „  jmsstrahlung  loss  parts  and  the  non- 
bremsstrahlung  part  of  the  wave  function,  and  clarifies  the  basic  role  of 
quantum  exchange  between  identical  particles  in  the  Quantum  1/f  Effect. 
This  present  formulation  is  briefly  presented  in  Sec.  11.2  and  11.3,  and  is 
described  in  more  detail  in  a  review  of  the  starting  points  of  the  Quantum 
1/f  Theory  submitted  to  Physical  Review.  Strictly  speaking,  the 
derivation  presented  in  Sec.  11.2  and  11.3  calculates  the  quantum  1/f 
fluctuations  in  space  along  the  scattered  beam  in  a  scattering  experiment, 
and  these  fluctuations  are  translated  into  fluctuations  in  time  by 
considering  the  passage  of  the  waves  past  a  fixed  point  in  space.  A  direct 
derivation  in  time  and  space  is  presented  in  Sec.  11.5  for  N-2,  and  yields 
the  same  result,  thereby  justifying  the  translation  of  wave  number 
spectra  into  frequency  spectra,  applied  in  Sec.  11.2-3. 
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Important  progress  was  achieved  through  an  extension  of  the  second 
-quantized  derivation  to  include  the  cross-correlations,  as  shown  in  Sec. 
11.6.  These  are  correlations  between  the  quantum  1/f  fluctuations  which  a 
physical  cross  section  or  process  rate  presents  for  scattering  in  a  certain 
direction,  and  the  similar  fluctuations  observed  on  the  same  cross  section 
or  process  rate  in  another  direction.  Knowledge  of  these  cross¬ 
correlations  is  necessary  for  the  calculation  of  the  resultant  mobility 
fluctuations  in  condensed  matter.  A  preliminary  calculation  of  the 
resulting  mobility  fluctuations  was  published  by  Kousik  et  al.9  in  1986, 
yielding  good  agreement  with  the  experiment,  and  is  now  being  revised  and 
refined  by  us  on  the  basis  of  these  more  rigorous  cross-correlation 
results. 

Another  important  problem  in  which  significant  progress  was  made, 
concerns  the  effect  of  a  finite  mean  free  path  on  quantum  1/f  noise  in 
condensed  matter.  In  principle,  the  Quantum  1/f  Effect  affecting  a  cross 
section  in  condensed  matter  should  differ  from  the  effect  calculated  for 
the  same  cross  section  in  isolated  conditions.  The  calculation  presented 
in  Sec.  11.7  shows,  indeed,  some  differences  expressed  in  terms  of  a 
correction  factor  (l-m2)*1  which  takes  into  account  multiple  scattering 
with  a  memory  coefficient  0<m<1  in  each  individual  scattering  process. 
However,  the  correction  factor  turns  out  to  be  1  in  the  case  of  ideally 
randomizing  collisions  (m»0)  which  erase  completely  the  memory  of  the 
incoming  particle  momentum  (e.g.,  lattice  scattering),  and  which  provide 
most  of  the  quantum  1/f  noise  anyway.  This  proves  that  all  our  earlier 
calculations,  which  were  equivalent  with  neglecting  this  correction,  did 
not  overestimate  the  quantum  1/f  noise,  and  were  an  excellent 
approximation,  in  spite  of  the  skepticism  of  some  critics. 

The  characteristic  functional  is  known  to  provide  the  most  complete 
description  of  a  random  process.  The  characteristic  functional  of  quantum 
1/f  noise  was  calculated  in  1982.  However,  the  Quantum  1/f  Effect  also 
affects  the  thermal  noise,  introducing  a  non-Gaussian  component.  Under 
the  present  Grant  the  author  was  able  to  calculate  the  characteristic 
functional  of  physical  thermal  noise  for  the  first  time,  as  we  briefly 
report  in  Sec.  11.8. 
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An  important  problem  related  to  1/f  nose  in  small  and  ultrasmall 
devices  is  the  transition  between  conventional  quantum  1/f  noise  which  is 
applicable  to  ultrasmall  devices,  and  coherent  state  quantum  1/f  noise 
which  is  applicable  in  the  limit  of  large  devices.  Due  to  the  complexities 
of  calculating  infrared  radiative  corrections  for  correlations  in  many- 
body  systems  in  condensed  matter  conditions,  the  progress  achieved  on 
this  problem  during  the  grant  period  was  unsatisfactory  and  limited  to  the 
creation  of  a  heuristic  interpolation  formula  based  on  the  author's 
physical  picture  which  associates  the  coherent  quantum  1/f  noise  with 
the  magnetic  field  of  the  collective  drift  motion  of  the  carriers.  While 
the  interpolation  formula  makes  sense  physically,  it  also  is  in  general 
agreement  with  the  experiment,  but  it  does  not  represent  a  mathematical 
solution.  This  interpolation  formula  is  presented  at  the  end  of  Sec.  11.9. 

The  whole  Quantum  1/f  Theory  can  be  reformulated  by  replacing 
photons  with  piezophonons  as  the  infraquanta.  This  also  includes  coherent 
state  quantum  1/f  noise.  Under  the  present  grant  the  coherent 
piezoelectric  quantum  1/f  noise  was  also  studied  and  is  briefly  presented 
in  Sec.  11.10.  More  experiments  are  needed  to  verify  it  in  detail,  although 
it  is  known  that  it  provides  in  principle  an  explanation  of  the  giant  Hooge 
parameters  observed  in  ferroelectric  materials. 

Along  with  the  further  development  of  the  quantum  1/f  theory,  a 
major  emphasis  was  placed  on  the  application  to  semiconductor  devices, 
which  is  presented  in  part  in  Sec.  III.  The  method  for  calculating  mobility 
fluctuation  Hooge  parameters  for  various  scattering  processes  in 
semiconductors  was  first  outlined  in  the  author's  1982  ONR  Final 
Technical  Report37.  This  method  was  applied  in  detail  by  Kousik  et  al.9 
and  resulted  in  good  agreement  with  the  experimental  data  for 
semiconductor  samples,  and  even  with  the  corresponding  Hooge 
parameters  estimated  without  averaging  over  the  energy  bands.  A  general 
introduction  to  the  application  to  pn  junctions  and  MIS  devices  is  provided 
in  Sec.  Ill .2. 

The  Quantum  1/f  Theory  was  first  applied  to  pn  junctions  by  van  der 
Ziel  and  Handel16  in  1985,  and  was  applied  to  n+p  junction  infrared 
detectors17  by  van  der  Ziel,  Handel,  Wu  and  Anderson  in  1986  and18  1989. 
This  later  work  is  presented  in  Sec.  III. 3,  and  led  to  the  formulation  of 
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optimal  design  principles  which  have  contributed  to  the  practical 
improvement  of  infrared  detectors.  In  1987  the  Quantum  1/f  Theory  was 
also  applied  by  the  author  to  MIS  infrared  detectors20,  as  shown  in  Sec. 

III.  4,  leading  to  analytical  formulae  which  allow  the  calculation  of 
quantum  1/f  noise  in  all  components  of  the  device  current  and  in  the 
resulting  detector  performance.  This  should  allow  for  practical  MIS 
device  optimization. 

In  cooperation  with  A.  van  der  Ziel  the  Quantum  1/f  Theory  was  used 
in  1986  to  calculate  the  Hooge  coefficient  of  Hgi-xCdxTe,  including  a 
relativistic  correction  in  the  Umklapp  scattering  part19.  The  correction 
is  needed,  because  the  very  low  effective  mass  of  electrons  in  narrow 
band  gap  Hgi-xCdxTe  leads  to  speeds  close  to  the  speed  of  light,  for  the 
crystal  momenta  usually  encountered  in  Umklapp  scattering  processes. 
Also  in  cooperation  with  van  der  Ziel,  the  quantum  partition  1/f  noise  in 
vacuum  pentodes21  was  reconsidered  and  shown  again  to  be  in  agreement 
with  the  experimental  data.  The  theory  was  also  successfully  applied  by 
van  der  Ziel  to,  and  verified  in,  secondary  emission22  and  photomultiplier 
tubes,  in  vacuum  photodiodes13,  silicon  pn  diodes23,  FETs13,  JfETs13 
BJTs24 

SQUIDs  are  important  high-tech  devices  limited  in  their  performance 
by  1/f  noise.  Therefore,  an  application  of  the  Quantum  1/f  Theory  to 
Josephson  junctions  and  SQUIDs  was  performed25,  based  on  the 
calculation  of  quantum  1/f  fluctuations  in  the  normal  resistance  value  of 
the  junction.  These  fluctuations  are  reduced  in  turn  to  the  Quantum  1/f 
Effect  present  in  the  elementary  cross  sections  which  define  the  normal 
resistance,  and  are  calculated  with  the  universal  quantum  1/f  formula. 

The  normal  resistance  enters  into  the  formula  of  the  critical  current  of 
the  Josephson  junctions.  Very  good  agreement  with  the  experiment  is 
obtained.  This  application  is  briefly  presented  in  Sec.  III. 5. 

The  few  selected  applications  presented  in  more  detail  in  this  report 
will  illustrate  the  practical  use  of  the  quantum  1/f  formulae,  and  the 
fruitful  interaction  between  theory  and  experiment  in  the  genesis  of  this 
new  discipline  linking  electrical  engineering  and  electrophysics  with 
quantum  electrodynamics. 
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II.  QUANTUM  1/f  NOISE  THEORY 

11.1.  SIMPLE  PHYSICAL  DERIVATION  OF  QUANTUM  1/f  NOISE 

The  Quantum  1/f  Effect,  also  known  as  conventional  quantum  1/f 
noise1*10,  is  a  fundamental  fluctuation  of  physical  cross  sections  and 
process  rates,  caused  by  the  infrared-divergent  coupling  of  current 
carriers  to  low  frequency  photons  and  other  infraquanta,  e.g.,  transversal 
phonons  with  piezoelectric  coupling,  or  electron  -  hole  pairs  on  the  Fermi 
surface  of  metals.  The  physical  origin  of  quantum  1/f  noise  is  easy  to 
understand.  Consider  for  example  Coulomb  scattering  of  current  carriers, 
e.g.,  electrons  on  a  center  of  force.  The  scattered  electrons  reaching  a 
detector  at  a  given  angle  away  from  the  direction  of  the  incident  beam  are 
described  by  DeBroglie  waves  of  a  frequency  corresponding  to  their 
energy.  However,  some  of  the  electrons  have  lost  energy  in  the  scattering 
process,  due  to  the  emission  of  Bremsstrahlung.  Therefore,  part  of  the 
outgoing  DeBroglie  waves  is  shifted  to  slightly  lower  frequencies.  When 
we  calculate  the  probability  density  in  the  scattered  beam,  we  obtain  also 
cross  terms,  linear  both  in  the  part  scattered  with  and  without 
bremsstrahlung.  These  cross  terms  oscillate  with  the  same  frequency  as 
the  frequency  of  the  emitted  bremsstrahlung  photons.  The  emission  of 
photons  at  all  frequencies  results  therefore  in  probability  density 
fluctuations  at  all  frequencies.  The  corresponding  current  density 
fluctuations  are  obtained  by  multiplying  the  probability  density 
fluctuations  by  the  velocity  of  the  scattered  current  carriers.  Finally, 
these  current  fluctuations  present  in  the  scattered  beam  will  be  noticed 
at  the  detector  as  low  frequency  current  fluctuations,  and  will  be 
interpreted  as  fundamental  cross  section  fluctuations  in  the  scattering 
cross  section  of  the  scatterer.  Although  the  wave  function  <p  of  each 
carrier  is  split  into  a  Bremsstrahlung  part  and  a  non-Bremsstrahlung  part, 
no  quantum  1/f  noise  can  be  observed  from  a  single  carrier.  A  single 
carrier  will  only  provide  a  pulse  in  the  detector.  Many  carriers  are  needed 
to  produce  the  quantum  1/f  noise  effect,  just  as  in  the  case  of  electron 
diffraction  patterns,  where  each  individual  particle  is  diffracted,  but 
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unless  we  repeat  the  experiment  many  times,  or  use  many  particles,  no 
diffraction  pattern  can  be  seen.  A  single  particle  only  yields  a  point  of 
impact  on  the  photographic  plate  in  diffraction,  or  a  pulse  in  the  detector 
in  the  1/f  noise  case.  While  incoming  carriers  may  have  been  Poisson 
distributed,  the  scattered  beam  will  exhibit  super-Poissonian  statistics, 
or  bunching,  due  to  this  new  effect  which  we  may  call  quantum  1/f  effect. 
The  quantum  1/f  effect  is  thus  a  many-body  or  collective  effect,  at  least 
a  two-particle  effect,  best  described  through  the  two-particle  wave 
function  and  two-particle  correlation  function. 

Let  us  estimate  the  magnitude  of  the  quantum  1/f  effect  by  starting 
with  the  classical  (Larmor)  formula  2q2a2/3c3  for  the  power  radiated  by 
a  particle  of  charge  q  and  acceleration  a.  The  acceleration  can  be 
approximated  by  a  delta  function  a(t)  =  Av8(t)  whose  Fourier  transform  Av 
is  constant  and  is  the  change  in  the  velocity  vector  of  the  particle  during 
the  almost  instantaneous  scattering  process.  The  one-sided  spectral 
density  of  the  emitted  Bremsstrahlung  power  2q2(Av)2/3c3  is  therefore 
also  constant.  The  number  2q2(Av)2/3hfc3  of  emitted  photons  per  unit 
frequency  interval  is  obtained  by  dividing  with  the  energy  hf  of  one 
photon.  The  probability  amplitude  of  photon  emission  [2q(Av)2/ 
3hfc3]1/2ei'r  is  given  by  the  square  root  of  this  photon  number  spectrum, 
including  also  a  phase  factor  e'T.  Let  \y  be  the  Schrodinger  wave  function 
of  the  scattered  outgoing  charged  particles,  which  can  be  constructed 
from  single-particle  wave  functions.  The  beat  term  in  the  probability 
density  Ivl2  is  linear  both  in  this  Bremsstrahlung  amplitude  and  in  the 
non-Bremsstrahlung  amplitude.  Its  spectral  density  will  therefore  be 
given  by  the  product  of  the  squared  probability  amplitude  of  photon 
emission  (proportional  to  1/f)  with  the  squared  non-Bremsstrahlung 
amplitude  which  is  independent  of  f.  The  resulting  spectral  density  of 
fractional  probability  density  fluctuations  is  obtained  by  dividing  with 
Ivl4  and  is  therefore 

lvl'4S|y|2(f)  »  8q2(Av)2/3hfNc3  =  2A/fN  -  r2Sj(f),  (1) 

where  a  *  e2/hc  =  1/137  is  the  fine  structure  constant  and  a  A  - 
4q2(Av)2/3hc3  is  known  as  the  infrared  exponent  in  quan+'im  field  theory, 
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and  is  known  as  the  quantum  1/f  noise  coefficient,  or  Hooge  constant,  in 
electrophysics. 

The  spectral  density  of  current  density  fluctuations  is  obtained  by 
multiplying  the  probability  density  fluctuation  spectrum  with  the  squared 
velocity  of  the  outgoing  particles.  When  we  calculate  the  spectral  density 
of  fractional  fluctuations  in  the  scattered  current  j,  the  outgoing  velocity 
simplifies,  and  therefore  Eq.  (1)  also  gives  the  spectrum  of  current 
fluctuations  Sj(f),  as  indicated  above.  The  quantum  1/f  noise  contribution 
of  each  carrier  is  independent,  and  therefore  the  quantum  1/f  noise  from  N 
carriers  is  N  times  larger;  however,  the  current  j  will  also  be  N  times 
larger,  and  therefore  in  Eq.  (1)  a  factor  N  was  included  in  the  denominator 
for  the  case  in  which  the  cross  section  fluctuation  is  observed  on  N 
carriers  simultaneously. 

The  fundamental  fluctuations  of  cross  sections  and  process  rates 
are  reflected  in  various  kinetic  coefficients  in  condensed  matter,  such  as 
the  mobility  p  and  the  diffusion  constant  D,  the  surface  and  bulk 
recombination  speeds  s,  and  recombination  times  x,  the  rate  of  tunneling  jt 
and  the  thermal  diffusivity  in  semiconductors.  Therefore,  the  spectral 
density  of  fractional  fluctuations  in  all  these  coefficients  is  given  also  by 
Eq.  (1).  This  is  true  in  spite  of  the  fact  that  each  carrier  will  undergo 
many  consecutive  scattering  processes  in  the  diffusion  process.  The 
quantum  1/f  noise  in  the  mobility  and  in  the  diffusion  coefficient  is 
practically  the  same8  as  the  quantum  1/f  noise  in  a  single  representative 
scattering  event  which  limits  the  mobility  or  the  diffusion  coefficient. 

Due  to  the  rapid  relaxation  of  concentration  fluctuations,  the 
quantum  1/f  fluctuations  of  scattering  cross  sections  will  only  be 
reflected  by  the  fluctuations  of  the  mobility  and  the  diffusion  constant  of 
the  carriers,  and  not  by  fluctuations  in  the  concentration  of  carriers. 

For  large  devices  the  concept  of  coherent  state  quantum  1/f  noise 
was  introduced11-  12.  In  this  case  the  Hooge  parameter  an  may  be  written 

«H  ■  («H)coh  *  2a/K  ■  4  6  10'3.  (2) 

where  a  «  1/(137)  is  the  fine  structure  constant.  This  is  of  the  same 
order  of  magnitude  as  the  empirical  value  an  -  2  10’3  that  Hooge  found  for 


long  devices.  It  is  therefore  proposed  that  Hooge's  empirical  value  for  an 
!s  due  to  coherent  state  quantum  1/f  noise,  so  that  it  has  a  very 
fundamental  origin. 

For  small  devices  (e.g.,  of  size  L  <  10  p)  we  apply  conventional,  or 
incoherent1*10,  quantum  1/f  noise  which  is  just  the  cross  section 
fluctuation  introduced  above  in  Eq.  (1).  In  that  case  an  may  be  written 

<*H  -  (aH)incoh  -  (4a/3*)[(Av)2/(C2)],  (3) 

where  Av  is  the  change  in  the  velocity  of  the  carriers  in  the  interaction 
process  considered.  This  expression  holds  for  any  1/f  noise  source 
describable  by  fluctuating  cross  sections.  Since  usually  (Av/c)2  «  1, 
except  for  carriers  with  a  very  small  effective  mass,  we  now  have  < 

3.1  10*3.  This  may  explain  the  low  values  of  (in  the  range  of  an  =  10*5 

•  10*9)  for  very  small  devices.  In  between  one  can  introduce  a  parameter 
s  -  f(L/L0)  where  L0  is  a  characteristic  size  and  write12. 

«H  -  (OH)incoh[1/(1  +  s)]  +  (aH)coh[s/(1  +  s)J,  (4) 

with  s  «  1  for  L/L0  «  1  and  s  »  1  for  L/Lq  »  1 .  This  describes  the 
transition  from  Eq.  (2)  to  Eq.  (3)  when  one  goes  to  devices  with  smaller 
and  smaller  sizes.  A  suggestion  for  the  calculation  of  s  was  presented  by 
the  author  at  the  Rome  1985  Conf.  on  Noise  in  Physical  Systems  and  1/f 
Noise,  s  -  2e2N'/m*c2  -  N'*5.5  •10*13cm,  where  N'  is  the  number  of 
carriers  per  unit  length  of  the  sample  or  device  in  the  direction  of  current 
flow,  and  5.5*1 0-^cm  is  twice  the  classical  radius  of  the  electron. 
According  to  this  rough  approximation12,  L 0  -  100p.  for  samples  with  a 
concentration  c  of  carriers  of  1015cm*3and  varies  proportional  to  c*1/2. 

When  we  apply  Eq.  (1)  to  a  certain  device,  we  first  need  to  find  out 
which  are  the  cross  sections  which  limit  the  current  through  the  device, 
and  then  we  have  to  determine  both  the  velocity  change  Av  of  the 
scattered  carriers  and  the  number  N  of  carriers  simultaneously  used  to 
test  each  of  these  cross  sections  or  rates.  Then  Eq.  (1)  provides  the 
spectral  density  of  quantum  1/f  cross  section  or  rate  fluctuations.  These 
spectral  densities  are  multiplied  by  the  squared  partial  derivative  of  the 
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current,  to  obtain  the  spectral  density  of  fractional  device  noise 
contributions  from  the  cross  sections  and  rates  considered.  After  doing 
this  with  all  cross  sections  and  process  rates,  we  add  the  results  and 
bring  the  fine  structure  constant  a  as  a  common  factor  in  front.  This 
yields  excellent  agreement  with  the  experiment13  in  a  large  variety  of 
devices  and  physical  systems.  This  general  principle  is  illustrated  on 
examples  of  practical  applications  in  Sec.  Ill,  but  here  in  Sec.  II  we  briefly 
first  describe  the  present  formulation  of  the  theory. 

11.2  DERIVATION  OF  THE  THEORY  FOR  N  BOSONS 

The  simplified  description  of  quantum  1/f  noise  was  presented 
above  in  the  elementary  terms  of  Schrodinger's  statistical  catalogue 
model,  without  using  second  quantization.  This  approach  is  natural  in 
view  of  the  close  connection  between  this  new  effect  and  diffraction 
which  is  usually  treated  without  second  quantization,  in  the  statistical 
catalogue  model  based  on  the  single-particle  solution  of  the  Schrddinger 
equation,  normalized  to  the  number  of  particles  N.  Just  as  the 
superposition  of  elementary  phase-shifted  waves  allows  for  the  simplest 
and  most  intuitive  description  of  diffraction  through  a  slit,  the 
description  of  quantum  1/f  noise  in  terms  of  interference  beats  between 
slightly  frequency-shifted  scattered  partial  waves  with  bremsstrahlung 
energy  losses  will  always  provide  the  simplest  and  most  elementary 
quantitative  derivation  of  the  quantum  1/f  effect,  easily  accessible  even 
at  the  undergraduate  level. 

Below  we  now  present  the  derivation  of  the  Quantum  1/f  Effect  in  a 
general  form  which  determines  the  scattered  current  j  from  the 
observation  of  a  sample  of  N  outgoing  particles.  The  minimal  outgoing 
sample  for  defining  particle-particle  correlations  in  the  scattered  wave 
consists  of  two  particles,  and  therefore  the  effect  can  calculated  for  the 
case  of  two  outgoing  particles.  Since  the  general  derivation  also  yields  a 
1/N  factor  for  bosons  and  a  factor  1/(N-1)  for  fermions,  and  since  the 
simplifying  restriction  to  N-2  has  given  rise  to  some  misinterpretations, 
a  presentation  of  the  general  case  of  N  bosons  or  N  fermions  will  be  of 
interest  for  us  as  a  generalization  of  the  results  derived  for  N-2  earlier 
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[#16,p.90]  on  the  pair  correlation  function.  In  this  Section  we  consider 
the  case  of  bosons. 

We  start  with  the  expression  of  the  Heisenberg  representation  state 
|S>  of  N  identical  bosons  of  mass  M  emerging  at  an  angle  0  from  some 
scattering  process  with  various  undetermined  bremsstrahlung  energy 
losses  reflected  in  their  one-particle  waves  <pj(£j) 

|S>  -  (N!)-l/2  Ili  JtPtwrfWir+Rl)  |0>  -  ni  Jd35i<pi(5i)  |S°>,  (5) 

where  y+(£j)  is  the  field  operator  creating  a  boson  with  position  vector  ^ 
and  |0>  is  the  vacuum  state,  while  jS°>  is  the  state  with  N  bosons  of 

position  vectors  £,•  with  i  -  1 . N.  All  products  and  sums  in  this  Section 

run  from  1  to  N,  unless  otherwise  stated. 

To  calculate  the  particle  density  autocorrelation  function  in  the 
outgoing  scattered  wave,  we  need  the  expectation  value  of  the  operator 

0(xi,x2)  -\y+(xi)\y+(x2)\}r(X2)v(xi),  (6) 

known  as  the  operator  of  the  pair  correlation.  This  operator  corresponds 
to  a  density  autocorrelation  function.  The  presence  of  two-particle 
coordinates  in  the  operator  O  does  not  mean  that  we  are  considering  two- 
particle  interactions,  it  only  means  that  the  expectation  value  which  we 
are  calculating  depends  on  the  relative  position  of  the  particles.  Using 
the  well  known  commutation  relations  for  boson  field  operators 

V(x)v+(y)  -  y+(y)y(x)  -  5(x  -  y),  (6'a) 

v(x)v(y)  -  v(y)v(x)  -  o,  (6'b) 

v+(x)y+(y)  -  v+(y)v+(x)  -  o,  (6*c) 

we  first  calculate  the  matrix  element: 


NkSP)qS°> 

*  S  jjlvS  mn8(i1v— Xi)5(T]j1-X2)5(^n~Xi)S(^m~X2)Z(j(j)n  ijS(rjj-^j),  (7) 
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where  |S°>  is  the  state  with  well  defined  particle  coordinates.  Here  the 
prime  excludes  p=v  and  m-n  in  the  summations  and  excludes  i-m,  i-n,  j-ji 
and  j=v  in  the  product.  The  summation  E(jj)  runs  over  all  permutations  of 
the  remaining  N-2  values  of  i  and  j.  On  the  basis  of  this  result  we  now 
calculate  the  complete  matrix  element 

<S|0|S>  -  [1  /N(N-1 )]  iVvS'mnfd^^TlvfdHmJdHn 

<Pn'(Tln)<Pv*(71v)<Pm(U)<Pn(^n)6(Tiv-Xi)8(iiu-X2)5(4n-Xi)5(4m-X2) 

■  [1/N(N-1)]  Z'jxvS'mn  <Pp.  (X2)<Pv  (Xi)(pm(Xi)(pn(X2).  (8) 

The  one-particle  states  are  spherical  waves  emerging  from  the 
scattering  center  located  at  x*0: 

<p(x)  *  (C/x)e'Kx  [i  +  2kib(k,l)e-iqx  a+k,i].  (9) 

Here  C  is  an  amplitude  factor,  K  the  boson  wave  vector  magnitude,  b(k,l) 
the  bremsstrahlung  amplitude  for  photons  of  wave  vector  k  and 
polarization  I,  while  a+k,i  is  the  corresponding  photon  creation  operator, 
allowing  the  emitted  photon  state  to  be  created  from  the  vacuum  if  Eq.  (9) 
is  inserted  into  Eq.  (8).  The  momentum  magnitude  loss  hq=Mck/Ks2jcMf/K 
is  necessary  for  energy  conservation  in  the  Bremsstrahlung  process. 
Substituting  Eq,  (9)  into  Eq.  (8),  we  obtain 

<S|0|S>  -  |C/x|4{N(N-1)  +  2(N-1)2k,i|b(k,l)|2[l+cosq(xi-x2)]},  (10) 

where  we  neglected  a  small  term  of  higher  order  in  b(k,l).  To  perform  the 
angular  part  of  the  summation  in  Eq.  (10),  we  calculate  the  current 
expectation  value  of  the  state  in  Eq.  (9),  and  compare  it  to  the  well  known 
cross  section  without  and  with  bremsstrahlung 


i  -  (hK/Mx2)[1  +  Iki|b(k,l)|2]  .  J0[1  +  aA/df/f], 


(11) 
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where  the  quantum  fluctuations  have  disappeared,  a=e2/hc  is  the  fine 
structure  constant,  aA«(2a/37t)(Av/c)2  is  the  fractional  bremsstrahlung 
rate  coefficient,  also  known  in  QED  as  the  infrared  exponent,  and  where 
the  1/f  dependence  of  the  bremsstrahlung  part  displays  the  well-known 
infrared  catastrophe,  i.e.,  the  emission  of  a  logarithmically  divergent 
number  of  photons  in  the  low  frequency  limit.  Here  Av  is  the  velocity 
change  h(K-K0)/M  of  the  scattered  boson,  and  f-ck/27c  the  photon 
frequency.  Although  the  bremsstrahlung  rate  described  by  Eq.  (11)  is 
generally  known,  we  derive  it  again  in  Sec.  11.4  below.  Eq.  (10)  thus  gives 

<S[0|S>  -  |C/x|4{N(N-1)  +  2(N-1)aAj[1+cosq(xi-x2)]df/f},  (12) 

which  is  the  pair  correlation  function,  or  density  autocorrelation  function 
along  the  scattered  beam  with  df/f=dq/q.  The  spatial  distribution 
fluctuations  along  the  scattered  beam  will  also  be  observed  as 
fluctuations  in  time  at  the  detector,  at  any  frequency  f.  According  to  the 
Wiener-Khintchine  theorem,  we  obtain  the  spectral  density  of  fractional 
scattered  particle  density  p,  (or  current  j,  or  cross  section  o)  fluctuations 
in  frequency  f  or  wave  number  q  by  dividing  the  coefficient  of  the  cosine 
by  the  constant  term  N(N-1): 

P"2Sp(f)  -  j'2Sj(f)  -  a-2SCT(f)  -  2aA/fN,  (13) 

where  N  is  the  number  of  particles  or  current  carriers  used  to  define  the 
current  j  whose  fluctuations  we  are  studying.  Quantum  1/f  noise  is  thus  a 
fundamental  1/N  effect.  The  exact  value  of  the  exponent  of  f  in  Eq.  (13) 
can  be  determined  by  including  the  contributions  from  all  real  and  virtual 
multiphoton  processes  of  any  order  (infrared  radiative  corrections),  and 
turns  out  to  be  aA-1,  rather  than  -1,  which  is  important  only 
philosophically,  since  aA«1.  The  spectral  integral  is  thus  convergent  at 
f-0. 

For  fermions  we  repeat  the  calculation  replacing  in  the  derivation  of 
Eq.  (10)  the  commutators  of  field  operators  by  anticommutators,  which,  as 
we  prove  in  detail  in  the  next  Sec.  11.3,  finally  yields  in  the  same  way 
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p-2Sp(f)  -  j-2Sj(f)  -  o-2Sa(f)  -  2aA/f(N-1 ),  (14) 

which  causes  no  difficulties,  since  N^2  for  particle  correlations  to  be 
defined,  and  which  is  practically  the  same  as  Eq.  (13),  since  usually  N»1. 
Eqs.  (13)  and  (14)  suggest  a  new  notion  of  physical  cross  sections  and 
process  rates  which  contain  1/f  noise,  and  express  a  fundamental  law  of 
physics,  important  in  most  high-technology  applications13. 

We  conclude  that  the  conventional  quantum  1/f  effect  can  be 
explained  in  terms  of  interference  beats  between  the  part  of  the  outgoing 
DeBroglie  waves  scattered  without  bremsstrahlung  energy  losses  above 
the  detection  limit  (given  in  turn  by  the  reciprocal  duration  T  of  the  1/f 
noise  measurement)  on  one  hand,  and  the  various  parts  scattered  with 
bremsstrahlung  energy  losses;  but  there  is  more  to  it  than  that:  exchange 
between  identical  particles  is  also  important.  This,  of  course,  is  just  one 
way  to  describe  the  reaction  of  the  emitted  bremsstrahlung  back  on  the 
scattered  current.  This  reaction  thus  reveals  itself  as  the  cause  of  the 
quantum  1/f  effect,  and  implies  that  the  effect  can  not  be  obtained  with 
the  independent  boson  model.  The  effect,  just  like  the  classical 
turbulence-generated  1/f  noise14,  is  a  result  of  the  scale-invariant 
nonlinearity  of  the  equations  of  motion  describing  the  coupled  system  of 
matter  and  field.  Ultimately,  therefore,  this  nonlinearity  is  the  source  of 
the  1/f  spectrum  in  both  the  classical  and  quantum  form  of  the  author's 
theory.  We  can  say  that  the  quantum  1/f  effect  is  an  infrared  divergence 
phenomenon,  this  divergence  being  the  result  of  the  same  nonlinearity. 

The  new  effect  is,  in  fact,  the  first  time-dependent  infrared  radiative 
correction.  Finally,  it  is  also  deterministic  in  the  sense  of  a  well 
determined  wave  function,  once  the  initial  phases  y  of  all  field  oscillators 
are  given.  In  quantum  mechanical  correspondence  with  its  classical 
turbulence  analog,  the  new  effect  is  therefore  a  quantum  manifestation  of 
classical  chaos  which  we  can  take  as  the  definition  of  a  certain  type  of 
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11.3  DERIVATION  FOR  THE  CASE  OF  N  FERMIONS 

In  the  case  of  fermions  the  calculation  is  similar,  except  for  the  use 
of  anticommutators  for  the  fermion  field  operators.  In  order  to  emphasize 
the  independence  of  our  results  on  the  representation  used,  and  to  show 
directly  how  the  calculations  presented  in  this  Section  can  be  performed 
without  second  quantization,  we  give  here  the  direct  calculation  in  terms 
of  a  Slater  determinant  for  the  state  of  N  scattered  fermions 


V . .  -  (N!)-i/2 


*i.(ri)  4>i.(r2)... 

4>i.(ri)  ♦b(r2)... 

...^(rn) 

‘Un)  <Ur2)... 

..‘Urn) 

(15) 


Here  r  combines  the  position  vector  x  and  the  spin  variable  s  for  each  of 
the  particles.  The  pair-correlation  function  is  obtained  by  integrating 
with  respect  to  the  coordinates  of  all  but  two  of  the  fermions 


A(ri ,r2)  -  /d3r3...d3r„<|¥i1...iN(ri...rn)|2>.  (16) 

Here  the  integrals  also  include  summations  over  the  spins;  the  expectation 
value  is  with  respect  to  the  phases  present  in  the  bremsstrahlung  parts  of 
the  wave  functions.  If  the  emitted  photons  are  included  in  the  final  state, 
the  expectation  value  is  also  taken  on  the  vacuum  background  of  the 
photons.  Assuming  orthonormality  of  the  functions  we  obtain 

N 

A(n,r2)  -  [1/N(N-1>]  £<|'fm(ri)'l'„(r2)-'Cm(r2)'t'n(ri)|2>.  (17) 

m;n»1 


To  display  the  spin  variables  explicitly,  we  write  <f>m(»’)  *  Xm(x)|s>  and  get 
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N/2 

ASS*(X1,X2)  -  [1/N(N-1 )]  £<[Xm*(Xl)<s|Xn*(X2)<s’|  - 

m;n-1 

Xm*(x2)<S,|Xn*(x1)<s|3  [Xm(x1)s>JXn{x2)s,>|  -  Xm(x2)s’>|Xn(x1)s>|]> 

N/2 

-  [1/N(N-1)1  X<[|Xm(xi)|2|X„(x2)|2  +  |Xm(x2)|2|Xn(xi)|2 

m;n-1 

-  Xm-(xi)Xm(x2)|<s|s'>|2X„-(x2)X„(x1)  -  (x,»x2)]>  (18) 


Here  the  symbol  (xi<=*X2)  designates  the  immediately  proceeding  term 
with  xi  and  X2  interchanged.  Considering  all  spin  orientations,  we  obtain 


N/2 

A(x,,x2)  -  An  +  At |  -  [1/N(N-1)J  2^  <[4lXml2|X„l2 

m;n-1 

-  Xm,(xi)Xn(x2)Xm'(x2)X„(xi)  -  (xi»x2)]>  (form  1) 


-[(|C|4)/xi2x22N(N-1)]{N2(1+Zk,||b(k,l)|2]2 

N/2 

-  2  X  exp[iKm(xi-x2)][l  +  2k,i|bm(k,l)l2exp[-iqm(xi-x2)] 

m-1 

N/2 

x  X  exp[-iKn(x1-x2)][1  +  SkMibn(kM,)|2exp[iqn(xi-X2)]}  (form  2) 
n-1 

-  t(C4)/x,  2X22(N-1  )]{N[1  +  Zk,l|b(k,l)|2]2 
N/2 

-  (2/N)  X  (I  +  aSk.HMk.lMScosqnfxnx^ 

n*1 

+  2k.i;kM,|bn(k,l)|2|b„(k,1l')l2eos(qn-q'n)(xi-x2)]J  (form  3) 


.((C4)/xi2X22(N-1)](N[1+Zk.l|b(k,l)|2]2-1-2Z|b(k,l)|2cosq(x,-x2)])  (19) 


This  form  of  the  pair-correlation  function  includes  the  1/N  factor  which 
multiplies  the  variable  (noise)  part.  The  crucial  point  in  the  derivation  of 
the  1/N  factor  was  an  elimination  of  the  rapidly  oscillating  terms 
expi(Kn-Km)(xi-X2)  with  Kn*Km  present  in  the  second  form  of  Eq.  (19) 
above,  an  elimination  indicated  through  the  approximation  sign  connecting 
the  second  form  to  the  third  form  above.  Indeed,  since  Km  differs  from  Kn 
by  much  more  then  the  momentum  change  corresponding  to  the  emission  of 
an  infraquantum,  these  terms  will  have  a  very  fast  oscillation,  and  will 
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not  yield  any  low  frequency  noise.  Since  they  are  also  small  in  magnitude, 
they  are  negligible.  This  provides  the  important  reduction  of  the  noise 
term  by  a  factor  N.  In  the  last  form  of  Eq.  (18)  we  note  that  the  first  two 
terms  are  constant  and  large,  and  do  not  yield  any  rapid  oscillations  which 
would  justify  elimination  of  any  cross  terms  with  Kn*Km.  Thus,  in  Eq.  (19) 
the  constant  part  of  the  pair  correlation  is  not  reduced  in  magnitude,  i.e., 
only  the  noise  part  is  affected. 

In  Eq.  (19)  we  have  used  again  the  form  of  the  single-particle  wave 
functions  given  by  Eq.  (9)  and  derived  below  in  Eq.  (41),  with  independent 
sets  of  phases  present  in  the  bremsstrahlung  energy  loss  parts  of  each 
particle.  Had  we  used  a  random  time-shift  in  each  of  the  single-particle 
wave  functions,  i.e.  just  a  random  initial  time  constant,  or,  equivalently,  a 
random  space-shift-constant,  all  results  would  have  been  exactly  the 
same.  In  fact,  the  random  shift  should  better  describe  the  initial  Poisson 
distribution  of  the  incoming  particles  which  are  scattered.  As  we  will 
see  at  the  end  of  Sec.  11.4  below,  the  random  phase  set  in  the 
bremsstrahlung  energy  loss  parts  is  the  set  of  random  initial  phases  of 
the  electromagnetic  field  oscillators.  Therefore  it  should  come  in  the 
same  way  for  all  particles.  However,  the  random  shift  will  eliminate 
expectations  of  the  cross  products  bm*(M)bn(k,l)  with  m*n  just  as  the 
sets  of  random  phases  used  by  us  did;  if  p  is  for  instance  a  random  space 
shift,  these  cross  terms  yield  contributions  to  the  pair-correlation 
function  of  the  form  <|b(k,l)|2  cosq(xi-X2  +  p)>  -  0,  where  the  average  is 
with  respect  to  p.  We  present  this  observation  here  as  an  afterthought, 
because  the  sets  of  random  phases  generated  by  a  shift  in  time  or  space 
will  appear  to  be  random,  but  will  still  contain  some  correlations.  The 
point  we  make  here  is  that  these  correlations  have  no  effect  on  our 
calculations,  so  we  can  continue  to  use  random  phases  for  our  purpose. 
There  may  be,  perhaps,  some  differences  in  higher-order  correlations 
which  we  do  not  consider  here. 

In  the  last  form  of  Eq.  (19)  we  have  neglected  the  higher-order  term. 
Performing  the  sum  in  the  same  way  as  in  Eqs.  (1 0)-(1 2),  or  as  we  show  in 
Sec. 11.4,  we  can  write  the  pair-correlation  function  for  fermions  in  the 
form 
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A(x-|  ,x2)  -[(|C|4)/xi2x22]{1  -2(N-1)*12k,i|b(k,l)|2[N-cosq(xi-x2)]}  (20) 

Dividing  again  the  variable  part  by  the  constant  term,  and  neglecting  small 
constant  terms,  we  finally  obtain  for  the  fractional  spectral  density  of 
the  fermion  current  and  cross  section  fluctuations 

Sj(k)dk/j2  ,  2aAdk/k(N-1)  -  Sj(f)df/]2  ,  2aAdf/f(N-1 1.  (21) 

11.4  DERIVATION  OF  THE  SINGLE-PARTICLE  WAVE  FUNCTIONS 

The  derivation  of  the  Schrbdinger  field  used  in  11.1  -  11.3  will  be 
performed  in  this  Section  with  the  help  of  the  Green's  function  method 
similar  to  the  method  used  in  a  earlier  calculation  by  Kroll  and  Watson15, 
extended  by  the  author  to  the  case  of  interaction  with  all  electromagnetic 
modes  of  the  universe7. 

It  is  most  convenient  to  describe  the  electromagnetic  field  in  terms 
of  plane  waves.  The  vector  potential  is  taken  in  the  radiation  gauge  as 

A(r,t)  -  2k,i(h2c/L2<ok)i/2Uk  |[ak>l(t)e'*<r  +  a‘k,,(t)e-*'].  (22) 

The  polarization  vectors  ukii  and  uk,2  are  mutually  orthogonal  unit 
vectors  perpendicular  to  k. 

The  Schrodinger  equation  for  an  electron  moving  in  a  vector 
potential  A  and  scattering  potential  V  is 

(1/2m)[-ibV  -  eA/cJ2\j/  +  V\y  -  itftjr.  (23) 

A  dot  has  been  used  to  indicate  the  time  derivative.  The  electromagnetic 
field  is  treated  as  a  classical  field  at  this  point.  In  order  to  eliminate 
the  A2  term  from  Eq.  (23),  we  write 

¥  -  exp[(-i/h)/t(e2/2mc2) A2dt']<t> .  (24) 

Thus,  Eq.  (4.2)  is  reduced  to 

[(-h2/2m)V2  +  (ieh/mc)A-V  +  V]0  -  ito<D. 


(25) 
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It  is  convenient  to  consider  first  the  influence  of  a  single  electromagnetic 
mode,  i.e.  a  single  term  from  Eq.  (22).  Therefore,  we  take  A  *  acos(cot+y), 
where  y  is  an  initial  phase  constant,  and  we  treat  VO  as  a  perturbation 
source  term.  The  solution  for  Eq.  (25)  is  an  incoming  plane  wave  plus 
scattered  waves,  given  by  the  integral  equation 

<J>ko  (r,t)  -  <t>ko  -  fd3xf  dt'GV(r')<Sko  (r'.f).  (26) 

Here  <|>ko  is  the  solution  of  the  homogeneous  equation,  i.e.  with  V  -  0,  and 
can  be  written  in  the  form 

Oko  -  exp[ik0r]exp[-(ih/2m)Jt(k20  -  2ek0*A/ch)dt].  (27) 

G  is  the  Green's  function  which  satisfies  the  equation 

[(-h2/2m)V2  +  (iehA-V/mc)  -  ih9/at]G  =  5(r-r')5(t-t').  (28) 

Given  A  *  acos(cot  +  y),  G  can  be  found  to  be 

G  *  [i/(2n)3b]Jd3keik(r’r,)exp{-(ih/2m)[k2t  -  2ek-asin(tot  +  y)/bcco]} 
x  exp[(ih/2m)(k2t’  -  2ek-asin(a>t  +  y)/bcco]  (29) 

In  the  first  Born  approximation  we  set  Ok  (r\t’)  =  <t>k  (r'.t*)  in  the  integral 
present  in  Eq.  (26)  and  obtain  for  the  scattered  wave 

Vs  -  [i/(2rc)3HjJd3x'/t  dt'V(r,)/d3keik(r-r,)exp{-(ih/2m) 

[k2t  -  2ek-asin(o)t  +  y)/hccu]}  x  exp[(ih/2m)(K2  -  ko2)t'] 
x  (exp[ie(k0  -  k)asin(o)t'  +  y)/2mcco]eikor,}•  (30) 

Using  the  relation 

eipsin(a)f+  y)  „  £  jn(p)ein{(ot'+Y)i  (31) 

n*-<«> 

where  Jn(p)  is  the  nth  order  Bessel  function,  we  expand  the  expression 
contained  in  curly  brackets  in  Fourier  series.  Then  Eq.  (26)  takes  the  form 
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<f>k  (r.t)  -  4>k  (r,t)  -  [-i^Tcja^JdSx'/tdt’Vfr^JdSkeifcfr-r  ) 

x  exp{-(ih/2m)[k2t  -  2ekasin(cot+Y)/hcco]} 
x  exp[(ifi/2m)(k2  -  k2)t’][  £  Jn(P)ein(«|,+T)]e,k  (32) 

n*-«o 

After  performing  the  integration  over  t‘  we  use  a  contour  integration 
method  for  k.  Then  Eq.  (32)  is  reduced  to 

«>k  (r.t)  -  <t>k  (r.t)  -  [-m/(2it)Ti2]  £  Jn(p)el"T[e»<(n)r/r] 

n*-<» 

exp{ih[k2(n)t-2ek(n)asin(cot  +Y)/hcco]/2m}  Jd3x,e'k(n)r,V(r,)eik'’,  (33) 

where 

p  »  -e(kn  -  k0)*a/mc£i)  -  -eQ-a/mchco,  (34) 

and  Q  is  the  momentum  transfer.  In  Eq.  (33)  k(n)  is  defined  by 

(hk(n))2/2m  -  (hk0)2/2m  -  nhco;  kn  *  k(n)  -  k(n)r/r.  (35) 

The  total  scattered  wave  can  be  written  as 

Vs  -  [-m/(2ji)h2r]  £  eik(n)rexp{-ih[k2(n)t-2ek(n)*asin(o)t+Y)/fico)]/2m} 

n»-«» 

X  Vk(n),k  Jn(P)e'nT  (36) 

where  Vk(n),k  -  Je-ik(n)-r’V(r')eiko*r,d3x'  is  the  scattering  matrix  element 
calculated  without  consideration  of  the  interaction  with  the 
electromagnetic  field  oscillators. 

So  far  the  electromagnetic  field  has  not  been  quantized  and  was 
considered  as  a  classical  field.  We  are  interested  in  the  corresponding 
expression  of  the  scattered  single-particle  wave  function  when  the 
electromagnetic  oscillators  are  quantized.  Therefore  we  first  linearize 
Eq.  (36)  with  respect  to  the  electromagnetic  potential  wherever  a 
dependence  on  \  is  present: 

Vs  -  [-m/(27t)h2r]VK,k  eiKr-'Et/hexp[ieK»asin(©t  +  Y)/cmco] 

{1  +  ei(wt-qr+7  )p/2  -  e-i(o>t-qr+y)|5/2} 

-  [-m/(27t)h2r]VK,k  ejKr-iEt/hexp[1  +  eiK-asin(cot  +  Y)/cmco] 

{1  +  e’(<ot-ci,'+'r)eQ*a/2mcho)  -  e-i(wt-<?r+y)eQ‘a/2mcfia)}. 


(37) 
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Since  k  -  (v/c)q«q,  the  r  dependence  of  the  elctromagnetic  potential  can 
be  neglected,  as  it  is  in  exp  [ieK*asin(cot  +7)/cma]. 

Here  we  have  introduced  the  notations  E  -  (toK)2/2m  and  K  *  k(0).  In  the 
last  form  we  will  neglect  the  term  with  sin(cot  +  y)  because  it  corresponds 
to  a  coherent  quantum  1/f  noise  contribution  which  has  been  considered 
elsewhere  before11-12,  and  because  this  term  can  be  considered  constant 
of  negligible  magnitude  in  the  calculation  of  equal  time  spatial 
correlations.  We  also  neglected  the  small  difference  between  k(1)  ■  K  - 
coK/vK  and  K  in  0  and  Vk(i),k  •  We  conclude  that  the  quantization  of  the 
electromagnetic  field  transforms  \ys  into  an  operator,  just  because  ak,i  is 
an  operator: 

<>(r,t)  -  (C/r)e'Kr-'Et/ti{i  -  Ik,|b*(k,l)expi(-a)t  +qr  -yOak.i 

+  Xk,|b(k,l)expi(tot  -  qr  +  Yi)a+k,i}-  (38) 

Here  we  have  introduced  b(k,l)  -  (1/2)0  and  the  constant  C  which 
designates  the  factor  in  front,  and  a  sum  which  includes  all 
electromagnetic  modes  with  annihilation  operators  ak.i-  In  the  space  of 
the  electron  states  <f»(r,t)  is  just  a  single-particle  wave  function.  We  have 
denoted  by  q  the  small  decrease  in  the  particle  momentum  required  by  Eq. 
(35).  We  have  q  =(3K/3E)hoj  -ck/v  »co/v,  with 

(L/2jt)3.4Tc.Z|<|b(k,l)|2>k2dk 

»  (e2Q2a2/4m2c2h2co2)*2k2dco/3c  -  aAdco/co,  (39) 

where  a  -  e2/4jthc  -  1/137,  A  -  (2Q2/37cm2c2),  and  where  <  >  is  an  angular 
average.  We  have  considered  the  spontaneous  emission  caused  by  vacuum 
fluctuations  only,  yielding 

<ak,l>  *  c(2h/coL3)1/2.  (40) 

The  annihilation  part  is  included  in  Eq.  (38),  but  does  not  contribute  to  the 
quantum  1/f  noise  on  the  background  of  the  electromagnetic  vacuum,  so 
we  can  simply  write 
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<p(x)  -  (C/r)e'Kr-iEt/h  [i  +  Zklb(k,l)e'M  -  qr  +  yi)  a+k (41) 

which  is  identical  with  Eq.  (9)  for  t=0,  which  in  turn  corresponds  to  the 
Heisenberg  representation  state  vector.  If  the  calculation  is  performed  on 
the  thermal  radiation  background,  however,  we  get  from  the  a+k,i  terms  a 
white  noise  contribution  added  to  the  quantum  1/f  noise;  the  latter, 
however,  remains  the  same4.  Here  we  have  performed  the  transition  from 
just  one  electromagnetic  mode  to  the  general  case  with  all 
electromagnetic  modes  adhoc,  but  in  our  paper7  this  transition  was 
presented  in  detail. 

In  Eq.  (38)  we  notice  the  presence  of  the  random  phases  y,  which 
where  introduced  as  initial  phases  of  the  electromagnetic  oscillators,  and 
which  are  independent  for  each  electromagnetic  mode  i  (of  given 
polarization  I  and  wave  vector  k)  of  the  universe.  Since  the  various 
scattered  particles  are  independent  of  each  other  and  of  the 
electromagnetic  modes,  we  have  to  consider  the  set  of  random  phases  y 
different  and  independent  in  the  wave  function  of  each  particle. 

11.5  DERIVATION  OF  THE  PAIR-CORRELATION  IN  TIME  AND  SPACE 

Working  in  the  Heisenberg  representation  as  above,  we  can 
generalize  the  derivation  presented  above  by  including  two  different 
times  in  the  operator  of  the  pair-correlation,  although  this  requires  more 
calculation.  To  simplify  the  integrals,  we  consider  again  the  state  of  two 
outgoing  fermions  very  far  from  the  place  where  they  have  independently 
suffered  the  same  interaction,  so  that  the  outgoing  spherical  waves  can  be 
approximated  by  plane  waves.  Starting  therefore  with  plane  waves 
similar  to  the  spherical  waves  used  in  Eqs.  (38)  and  (9),  with  q  *  qK/K  = 
ckK/Kv  and  q'  =  ck'K/Kv,  we  obtain  again  the  pair-correlation  function  in 
the  form 

A(x-| ,  ti;  x2,  t2)  -  l/alda^T^'dan’ 

(exp(-iKiV)  +  Zk.iS-r(k,,s’)exp[-i(K-q’)Ti,]ak-,s1} 

{exp(-iK^)  +  Zktsb*(k,s)exp[-i(K-q)^]ak,s} 
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{exp(iK4)  +  Ik,sb(k,s)exp[i(K-q)51a+k,s} 

{exp(iKri)  +  Zk^s'P(k^s,)exp(i(K-q,)1^]a+k^s•} 

(1/4)Sss-<S0ssi0n  +  Oti  +  Oxt  +  Oxxlxt.ti;  Xz,h  |S°ss'  >•  (42) 

The  operator  of  the  pair-correlation  function  now  contains  two 
consecutive  times.  The  creation  and  annihilation  operators  for  particles 
obey  anticommutation  relations  similar  to  the  commutation  relations  in 
Eqs.  (6’a)-(6'c).  Using  these  operators,  the  field  operators  can  be  expanded 
in  terms  of  plane  waves, 

Vs(r.t)  -  V-i/2£p  expi(pr-o)t)CpS, 

V+s(r,t)  -  V-i/2£p  expi(-pr+o)t)c+pS  (43) 

which  also  contain  their  time  dependence,  as  needed  for  the  operator  of 
the  pair  correlation  function  in  the  Heisenberg  representation.  This 
operator  has  an  expectation  value  which  can  be  written,  for  spin  up  only, 
in  the  form 

<S°  |Ott(xi,  ti,  X2,  t2)|S°  > 

-  ,ti )v+(X2,t2)\|/(X2,t2)v(x1  ,ti )y+(S)V+(h)  |0> 

■  (1/V4)  £  5^  expi[m,T|,  “jx'X'i  +a)(p.')ti  -v'X2 +a)(v')t2 

mnm'n’nviiV 

+11X2  -«(n)t2  +vxi  -<B(v)ti  -mr\  -n^]<0|cm-cn-c+v'C+p’CpCvc+mc+n|0  > 

-  (M/h)6tr3t2-3exp{(iM/2h)[(xi-^)2/t1  +(x2-Ti)2/t2 

-(xi-n')2/tl  -(x2-5’)2/t2]} 

+(M/h)6ti  *3t2'3exp{(iM/2h)[(xi  -t]  )2/tt  +(x2-^)2/t2 
-(Xi-^)2/tl  -(X2-TT)2/t2]} 

-(M/h)6ti-3t2-3exp{(iM/2h)[(xi-^)2/t1  +(x2-n)2/t2 
-(Xi-V)2/tl  ~(X  2 “TJ  *) 2/t2 ] } 

-(M/h)6tr3t2-3exp{(iM/2h)[(xi-Ti)2/t1  +(x2-^)2/t2 
-(xi-Ti’)2/ti  -(X2-V)2/t2)}.  (44) 

Substituting  into  Eq.  (42),  and  integrating  with  respect  to  £ ,r\ and  r|',  we 
obtain  for  the  part  with  spin  up  only 
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8Att(xi,  ti;x2,  t2) 

-  <0{1  +2k,P*(k')expi[q,xi  -bq'(K  -q'^t-j/Mjak1} 

{1  +Ikb*(k)expi[qx2  -hq(K  -q/2)t2/M]ak} 

{1  +Skb(k)expi[qxi  +fcq(K  -q/2)ti/M]a+k} 

{1  +Sk'P(k')expi[q'x2  +hq'(K  -q72)t2/M]a+k*}0> 

-  1  +Sk’|b(k,)|2expi[q'(xi  -x2)  -bq'(K  -q'/2)(ti  -t2)/M] 

+Zk|b(k)|2expi[q(x2  -xi)  -fcq(K  -q/2)(t2 
+Zk,k'|b(k)b(k‘)!2expi[(q  -q’)(x2  -xi)  -hK(q  -q')(t2  -ti)/M 
+h(q2  -q'2)(t2  -ti)/2M],  (45) 

In  a  similar  way  the  expectation  values  corresponding  to  the  uspin 
orientations  and  the  two  identical  terms  with  all  spins  reversed  are 
calculated  and  integrated  with  respect  to  and  t|\  Putting  together 
the  terms  for  all  spin  orientations  as  before,  we  finally  obtain  the  desired 
pair  correlation  function  in  the  form 

A(x,,ti;x2,t2)  -  1/2  +  £it,i|b(k,l)|2{2  -  cosq[xi-x2  -v'(ti-t2)]) 
+£kkT|b(k,l)|2|b(k\l)|2{1  -(1/2)cos(q-q’)[x,-x2  -w"(tj-t2)]>;  (46) 

v'»  v(1  -q/2K)  -  v(1  -e/4E)  =»  v;  vM-  v[1  -(q+q')/2K]  -  v[1  -(e+e')/4E]  =  v. 

The  approximations  v*  ~  v  and  v"  =■  v  are  justified  because  the  soft  photon 
energy  e  -  4  10-15eV  for  1Hz  is  negligible  compared  to  the  energy  E  of  the 
particles  which  may  be  of  the  order  of  leV.  This  result  directly  proves 
the  validity  of  our  translation  of  wave-number  spectra  into  frequency 
spectra  in  Sec.  11.2  and  11.3.  The  extension  to  the  case  of  bosons  is  trivial. 

11.6.  QUANTUM  1/f  CROSS-CORRELATIONS  AND  SPECTRA 

The  calculation  of  quantum  1/f  noise  in  various  kinetical  coefficients 
governing  transport  in  condensed  matter,  requires  a  knowledge  of  the 
cross-correlation  of  the  outgoing  current  densities  scattered  by  any 
process  into  different  directions  with  wave  vectors  K’  and  K"  when  the 
incoming  particles  had  wave  vectors  Ki  and  K2.  The  particles  are 
assumed  to  be  identical,  of  mass  M,  with  wave  functions  which  overlap 
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somewhere.  As  we  have  seen  in  the  preceding  Secs.  11.2-3  on  the  general 
derivation  of  quantum  1/f  noise  in  physical  cross  sections,  the  time 
correlations  and  frequency  spectra  can  be  readily  obtained  from  the 
corresponding  spatial  correlations  along  the  outgoing  beams  in  the 
direction  of  the  detector  which  is  now  replaced  by  two  detectors  in  the 
directions  of  K'  and  K"  respectively.  This  is  what  we  called  translation 
from  wave-numbers  to  frequencies  at  the  end  of  the  last  Section.  Working 
again  in  the  Heisenberg  representation,  we  limit  ourselves  to  the  simplest 
case  of  two  particles  described  by  the  outgoing  spherical  waves 

<p(x)  -  (C/x)exp(iKi x)[1  +  Ik,|b(k,l)exp(-iq1x)a+k,i]10>  (47) 

x(x)  =  (C/x)exp(iK2x)[1  +  Ik,lP(k,l)exp(-iq2x)a+k,i]|0>  (48) 

where  now  b  corresponds  to  the  momentum  change  K'-Ki,  while  p  has  an 
independent  phase  and  corresponds  to  the  momentum  change  K"-K2  of  the 
second  particle.  As  in  the  preceding  Sec.  11.2,  for  bosons  we  obtain  the 
current  density  autocorrelation  function 

A(K1,K2,K'iK";xix2)  =■  <S|0|S>  =  (K'K72M)<0|x(xi)<p(x2)|2 

+5c'(Xr)q>‘(*2)x(*2)<P(Xl)  +X*CX2)<p*(xi)x(xi)<P(X2)  +lx(*2)<f>(*l)l2|0> 

=  (K'K7M){[1  +  Xk,l|b(k,l)|2][1  +  £w|p(k,l)|2] 

+(1/2)exp[iK,(x2-xi)-iK"(X2-x1)][1  +  Ek,llP(k,l)l2exp-iq(X2-xi)] 

[1  +  Ek.l|b(k,l)|2expiq(x2-xi)]  +  [xi  <->  x2J} 

-  (K'K7M){[1  +Ik.i|b(k,l)|2][l  +Zk,llP(k.l)|2]  +  cos(K'-K')(X2-xi) 

+  £k.llP(k,l)|2cos[(K,-K--q)(x2-x1)  +Ik,i|b(k,l)|2cos[(K'-K"+q)lx2-x1) 

+  Ik,i£kj'IP(k,l)|2|b(kM’)|2cos[(K'-K-+q'-q)(X2-xi)]}.  (49) 
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Except  for  K'=K"  all  terms  which  are  not  constant  are  rapidly  oscillating 
and  do  not  yield  low  frequency  noise  contributions;  as  usual,  only 
practically  unobservable  high  1/Af  noise  terms  are  present  in  this  case  in 
the  quantum  fluctuations.  For  K'»K",  however,  we  get  low  frequency  noise 
down  to  arbitrarily  low  frequencies,  because  the  momentum  magnitude 
changes  q=kMc/hK  and  q'=k'Mc/hK  caused  by  the  soft  photon  emission  are 
arbitrarily  small.  We  obtain  the  spatial  quantum  1/f  current  density 
cross-correlation 

A(Ki,K2,K\K-;xiX2)  -  (h2K’K7M){[1+2icii|b(k,l)|2][1+Ek,ilP(k,l)|2] 

+  {1  +  Xk.ll|b(k,l)|2  +  ||3(k,l)|2]cosq(X2-xi)}8K',K' 

+  2k,|Xk',llP(k,l)l2|b(kM')l2cos[(q'-q)(X2-x,)]5K',K-}-  (50) 


Note  that  the  restriction  indicated  by  the  Kronecker  symbol  refers  only  to 
the  magnitude  of  the  outgoing  wave  vectors,  or  to  the  kinetic  energies  of 
the  outgoing  particles,  and  that  their  outgoing  and  incoming  directions  can 
be  different.  Therefore,  only  particles  of  the  same  energy  have  quantum 
1/f  noise  cross-correlations,  and  particle  groups  of  different  energy  will 
yield  independent  quantum  1/f  noise  contributions,  a  result  which  was 
anticipated  empirically  by  Kleinpenning  in  semiconductors.  Performing 
the  angular  integrations  contained  in  the  softaron  lattice  sums  of  Eq.  (50), 
as  indicated  in  Eq.  (39),  and  using  the  translation  relation  dk/k  =  df/f, 
with  f*ck/2ic,  we  obtain 

A(K1,K2,K\K";xix2)  =  (h2K'K7M){[1  +aAiJdf/f][1  +<xA2Idf/f] 

+  [1  +  aJ(Ai  +  A2)cos  q(x2-xi)  (df/f)]5K\K" 

+  a2Ai  A2J(df/f)J(dfVf')cos[(q,-q)(x2-xi)]5K,,K"}-  (51) 

Here  Ai  -  2h2(K’-Ki)2/37cM2c2  and  A2  -  2fc2(K"-K2)2/3*M2c2.  We  may  call 
this  approximation  the  independent  bremsstrahlung  model.  Neglecting  the 
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last  term  which  is  a  noise  of  noise  contribution,  applying  the  Wiener- 
Khintchine  theorem,  and  using  again  the  relation  dk/k  *  df/f,  to  connect 
wave  number  spectra  to  frequency  spectra,  we  obtain  the  cross-spectral 
densities  of  the  rate  fluctuations  Aw(Ki.K')  from  Ki  to  K',  and  Aw(K2,K") 
from  K2  to  K" 

Saw(Ki,K2,K’,K")  -(a/2f)(Ai+A2)<w(Ki,K,)><w(K2,K-)>5K-,K-(2/N)  (52) 

for  bosons,  while  for  fermions  we  obtain  in  the  same  way 

SAw(Ki,K2,K,,K-)*(ayf)(A1+A2)<w(Ki,K’)><w(K2,K-)>5K-,K-(1/N-1).  (53) 

The  last  factors  have  been  added  in  an  attempt  to  generalize  the  results, 
derived  here  for  N=2,  to  the  case  of  N  particles  in  the  final  state,  on  the 
basis  of  the  result  of  the  preceding  Sections  11.2  and  11.3. 

11.7  EFFECT  OF  A  FINITE  MEAN  FREE  PATH  ON  QUANTUM 

1/f  NOISE 


The  very  low  frequency  radiation  emitted  by  current  carriers  in 
condensed  matter  is  important  for  the  calculation  of  quantum  1/f  noise. 
Here  the  spontaneously  emitted  radiation  is  calculated  for  a  carrier 
subject  to  a  large  number  of  closely  consecutive  scattering  processes. 

The  result  of  the  calculation  is  similar  to  the  radiaiton  from  a  single 
representative  scattering  process,  with  a  correction  factor  of  the  order  of 
unity. 

In  this  Section  we  show  that  conventional  quantum  1/f  noise  from  a 
large  number  of  closely  consecutive,  random,  scattering  events  can  be 
calculated  from  the  average  properties  of  a  single  scattering  process,  and 
is  close  to  the  quantum  1/f  noise  of  an  isolated  scattering  cross  section 
in  the  case  of  randomizing  collisions. 

Here  we  calculate  the  extreme  low  frequency  (e.I.f.)  radiation  from  a 
large  number  N  of  very  frequent  consecutive  scattering  events  and 
compare  with  what  we  would  expect  from  a  single  isolated  scattering  of 
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the  current  carrier.  It  is  well  known  that  these  radiation  energy  losses 
also  determine  the  conventional  quantum  1/f  noise1*8. 

We  start  with  a  review  of  the  derivation  of  the  power  radiated  in  a 
collision  with  acceleration  v  of  the  carrier  of  charge  e: 


A  -  (~4jc/c)j,  (54) 

A(x,t)  -  (1/c)f{J[x\t  -  (|x  -  x ' |/c)]/| x  -  x'|}d3x' 

-  (1/cr)/j[(x\t  -jx*  -  x|/c)]d3x';  (55) 

E  -  -(1/c)A  -  ev/c2r,  B  -  (k/jk|)xE;  (56) 

Y  -  (c/4n)ExH,  P  -  /YdS  -  (2e2/3c3)v2.  (57) 


For  one  collision  we  separate  a  totally  random  part  of  zero  average  Av, 
and  an  average  part  proportional  with  the  previous  velocity,  in  the  total 
velocity  change  (Av)t0t: 

v  -  8(t)(Av),ot 

(Av)tot  -  Av  -  yv  (Definition  of  y);  0  <  y  <  1.  (58) 

We  note  that  for  randomizing  collisions  with  little  or  no  memory  the 
randomness  parameter  y  approaches  the  unity  from  below,  and  is  very 
close  to  unity.  For  N  consecutive  sudden  accelerations  we  obtain  for  v  a 
sum  of  N  delta  functions.  Expanding  the  delta  functions  as  Fourier 
integrals,  the  total  energy  emitted  by  one  particle  is 

W  -  J  Pdt  -  (2e2/3c3J  dt|[(Ai  v)totJ  e2,tiftdf 

+(A2V)t0J  e2*if(t-t  )df  +  . +(ANv)toJe2*if(t-t  >  ]|2.  (59) 

Neglecting  the  short  times  between  collisions,  i.e.,  setting  ti  -  t2  » . * 

tN-i  ■  0,  we  obtain  in  the  low  frequency  limit,  i.e.,  if  we  do  not  care  about 
errors  in  the  integrands  and  in  the  resulting  power  spectral  density  as 
long  as  they  are  restricted  to  the  high  frequency  region, 

W  -  (2e2/3c3)J  dtf  dfj  dfe2*'<M*)t 

|Aiv  *yv  +A2V  -y(Aiv  +(1-y)v]+A3V  -y{A2V  +(1-y)[Aiv  +(1-y)v]} 

+A4V  -y{A3V  +(1-y)[A2V  +(1-y)(A-|V  +(1-y)v)]}  +  . |2  forf-»0.  (60) 
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Summing  the  geometric  series  repeatedly, 


N-1  N-2 

W  -  (2e2/3c3)J  df|-yv  £(1-Y)n  +Aiv  [1  -y  Id-7)"] 

n-0  n-0 

N-3 

+  A2v[1  -  Y  X(1-Y)n]  + . I2 

n-0 

-  (2e2/3c3)J  df|v[1  -(1-y)n]  +Aiv(1-y)n'1  +A2v(1-y)n*2  +....+Anv|2.  (61) 

Finally,  we  take  the  square  and  we  average  over  the  independent 
random  Av  vectors.  Summing  up  geometric  series  again,  we  obtain  for  N 
going  to  infinity  in  the  limit 

<W>  -  (2e2/3c3)j  df{<v2>[1-(1-Y)N]2  +<(Av)2>[1  -(1-y)2N]/[1  -(1-y)2]} 

-»  (4e2/3c3)J  df{<v2>  +  <(Av)2>/y(2  -y)}  (62) 

The  e.I.f.  emission  of  radiation  from  a  large  number  of  statistically 
similar,  closely  spaced  consecutive,  random  collisions  of  a  current  carrier 
is  thus  close  to  the  emission  from  just  one  average  collision  or  scattering 
process,  as  I  had  implied  before. 

The  corresponding  quantum  1/f  effect  in  the  current  carried  just  by 
this  carrier  will  therefore  be1*9 

S8j(f)/(j2)  -  2(4e2/3c3hf){<v2>  +  <(Av)2>/y(2  -y)} 

-  (4a/37tf){<v2>/c2  +  <(Av)2>/c2y(2  -y)}.  (63) 

The  two  terms  in  curly  brackets  should  be  roughly  equal  in  average,  in 
thermal  equilibrium  conditions.  Indeed,  by  equating  the  average  speed  of 
the  particle  before  and  after  a  collision  <(v  -yv  +av)2>  -  <v2>,  we  obtain 
<v2>  »  <(Av)2>/[y(2  -  y)]. 

The  randomness  parameter  y  can  be  replaced  by  a  memory  parameter 
m  -  1  -  y,  where  again  0  <  m  <  1 .  Then  the  corrective  factor  becomes 
1/[y(2  -  y)]  -  1/[1  -  m2].  This  shows  that  the  corrective  factor  will  not  be 
less  than  unity. 
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For  condensed  matter  applications,  based  on  the  classical  concept  of 
many  consecutive  random  collisions  of  the  same  particle,  one  might 
question  the  low  frequency  content  of  the  delta  functions  used  above.  In  a 
quantum  discussion,  however,  we  can  not  identify  the  particles  in  general, 
and  all  we  can  say  is  that  N  scattering  processes  have  taken  place.  Due  to 
quantum  exchange  processes,  the  particles  may  have  been  an  indefinitely 
long  time  in  the  initial  state,  and  may  persist  in  the  final  state  an 
indefinite  time.  Even  in  the  classical  case,  however,  the  low  -  frequency 
content  of  the  delta  functions  used  in  this  Section  is  limited  only  by  the 
reciprocal  observation  time  "M,  because  the  velocity  of  the  carrier  is 
well  defined  over  the  whole  duration  T  of  the  experiment. 

Note  that  for  randomizing  collisions  with  little  or  no  memory  the 
corrective  factor  1/[y(2  -  y)]  is  close  to  unity.  This  proves  that 
conventional  quantum  1/f  noise  can  in  many  cases  be  applied  without 
including  the  finite  mean  free  path  correction. 

11.8  CHARACTERISTIC  FUNCTIONAL  OF  PHYSICAL  THERMAL  NOISE 
WHICH  INCLUDES  EQUILIBRIUM  1/F  NOISE 

In  a  published  paper$  the  characteristic  functional  of  quantum  1/f 
noise  was  derived.  This  result  was  applicable  to  quantum  1/f  noise  in  any 
cross  section  or  process  rate,  and  in  currents  or  voltages.  In  an  early 
paper2^  the  application  of  quantum  1/f  noise  to  thermal  equilibrium  noise 
was  performed,  and  the  quantum  1/f  Nyquist  theorem  was  formulated. 
Deviations  from  the  Gaussian  form  of  thermal  noise  were  expressed27  in 
terms  of  the  characteristic  function. 

The  present  Section  reports  the  final  step  in  the  study  of  thermal 
noise  with  infrared  radiative  corrections,  by  presenting  the  characteristic 
functional  of  physical  thermal  noise.  We  call  it  "physical"  because  we 
refer  to  the  actually  observed  thermal  equilibrium  (Johnson)  noise,  rather 
than  to  the  strictly  Gaussian  noise  expected  from  the  Nyquist  formula 
without  infrared  radiative  corrections.  The  small  quantum  1/f 
contributions  come  in  as  time-dependent  infrared  radiative  corrections, 
required  by  the  interaction  of  the  current  carriers  with  the 
electromagnetic  field,  or  by  the  reaction  of  the  bremsstrahlung  back  on 
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the  current  which  has  produced  it.  Knowing  the  characteristic  functional 
of  a  process  is  the  best  we  can  do,  the  highest  level  and  most 
comprehensive  and  complete  description  possible  for  a  random  process. 

In  terms  of  current  fluctuations  in  equilibrium,  to  derive  the 
characteristic  functional  of  the  physical  thermal  noise  variable  £,  we 
express  it  in  terms  of  the  unmodulated  theoretical  Nyquist  noise  variable 
x: 


\  -  x(G/G0)1/2  -  x  +  xy,  (64) 

where  y  ■  8G/G  is  the  fractional  fluctuation  in  the  conductivity  G  of  the 
conductor  whose  thermal  equilibrium  current  fluctuations  we  are 
considering.  Let  the  quantum  1/f  noise  variable  y  obey  a  Gaussian 
amplitude  distribution  of  dispersion  02  and  x  one  of  dispersion  01.  In 
terms  of  the  zero-order  Bessel  function  of  imaginary  argument,  the 
product  z=xy  will  have  an  amplitude  distribution  P(z)=(1/7caia2)K0(z/aia2) 
which  has  an  elementary  characteristic  function 

X(v)  -  (1  +  oi2o22v2)'1/2.  (65) 

Although  z  and  x  are  not  independent,  we  have  been  able  to  calculate  the 
characteristic  function  of  %  as  an  elementary  function  in  the  form 

X(k)  -  (1  +  k2ai2a22)'1/2  exp[-k20l2/2(1  +  k2ai2a22)],  (66) 

while  the  amplitude  distribution  itself  can  not  be  expressed  in  elementary 
functions. 

In  the  second  part  of  our  work  [#31,  p.91]  we  derive  the  characteris¬ 
tic  functional  of  quantum  1/f  noise,  and  we  familiarize  the  reader  with 
characteristic  functionals  in  general.  Finally,  using  also  Eq.  (66),  we 
derive  the  characteristic  functional  of  physical  thermal  noise  in  the  form 


F[k(f)]  -  exp{-(1/2)Jln[1  +  8kBTG(<xA/fN)(f/f0)aAk2(f)]df} 
exp{-2kBT  G  Jk2(f)[i  +  8kBTG(aA/fN)(f/f0)aAk2(f)]-idf}. 


(67) 
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Here  aA  is  the  infrared  exponent  <xA-(4a/3rc)<(Av/c)2>,  which  is  known  to 
enter  the  expression  of  both  the  basic  quantum  1/f  noise  formula  and  the 
characteristic  functional  of  quantum  1/f  noise  in  two  ways:  as  a 
coefficient  and  as  an  additional  exponent  of  the  frequency  f. 

Sommerfeld's  fine  structure  constant  a-e2/hc  -1/137  is  well  known,  and 
so  is  the  Boltzmann  constant  ke-  Eq.  (4)  gives  the  characteristic 
functional  of  physical  thermal  noise  for  a  sample  containing  N  current 
carriers  with  an  average  velocity  change  <(Av)2>  in  the  scattering 
processes  which  determine  their  mobility. 


11.9  COHERENT  STATES  AND  CONVENTIONAL  QUANTUM  1/f  EFFECT 

An  electrically  charged  particle  includes  the  bare  particle  and  its 
field.  The  field  has  been  shown  in  the  last  two  decades  to  be  in  a  coherent 
state,  which  is  not  an  eigenstate  of  the  Hamiltonian.  Consequently,  the 
physical  particle  is  not  described  by  an  energy  eigenstate,  and  is 
therefore  not  in  a  stationary  state.  In  this  Section  we  show  that  the 
fluctuations  arising  from  this  non  -  stationarity  have  a  1/f  spectral 
density  and  affect  the  ordered,  collective,  or  translational  motion  of  the 
current  carriers.  This  "coherent"  quantum  1/f  noise  should  be  present 
along  with  the  familiar  quantum  1/f  effect  of  elementary  cross  sections 
and  process  rates  introduced  ten  years  ago,  just  as  the  magnetic  energy  of 
a  biased  semiconductor  sample  coexists  with  the  kinetic  energy  of  the 
individual,  randomly  moving,  current  carriers.  The  amplitude  of  the 
quantum  1/f  effect  is  always  the  difference  of  the  coherent  quantum  1/f 
noise  amplitudes  in  the  "out"  and  "in"  states  of  the  process  under 
consideration  and  dominates  in  small  samples,  while  large  samples  should 
exhibit  the  larger  coherent  quantum  1/f  noise. 

A  physical,  electrically  charged,  particle  should  be  described  in 
terms  of  coherent  states  of  the  electromagnetic  field,  rather  than  in 
terms  of  an  eigenstate  of  the  Hamiltonian.  This  is  the  conclusion  obtained 
from  calculations28  of  the  infrared  radiative  corrections  to  any  process 
performed  both  in  Fock  space  (where  the  energy  eigenstates  are  taken  as 
the  basis,  and  the  particle  is  considered  to  have  a  well  defined  energy)  and 
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in  the  basis  of  coherent  states.  Indeed,  all  infrared  divergences  drop  out 
already  in  the  calculation  of  the  matrix  element  of  the  process 
considered,  as  it  should  be  according  to  the  postulates  of  quantum 
mechanics,  whereas  in  the  Fock  space  calculation  they  drop  out  only  a 
posteriori,  in  the  calculation  of  the  corresponding  cross  section,  or 
process  rate.  From  a  more  fundamental  mathematical  point  of  view,  both 
the  description  of  charged  particles  in  terms  of  coherent  states  of  the 
field,  and  the  undetermined  energy,  are  the  consequence  of  the  infinite 
range  of  the  Coulomb  potential29.  Both  the  amplitude  and  the  phase  of  the 
physical  particle's  electromagnetic  field  are  well  defined,  but  the  energy, 
i.e.  the  number  of  photons  associated  with  this  field,  is  not  well  defined. 
The  indefinite  energy  is  required  by  Heisenberg's  uncertainty  relations, 
because  the  coherent  states  are  eigenstates  of  the  annihilation  operators, 
and  these  do  not  commute  with  the  Hamiltonian. 

A  state  which  is  not  an  eigenstate  of  the  Hamiltonian  is 
nonstationary.  This  means  that  we  should  expect  fluctuations  in  addition 
to  the  (Poissonian)  shot  noise  to  be  present.  What  kind  of  fluctuations  are 
these?  This  question  was  answered  in  a  published  paper11.  The  additional 
fluctuations  were  identified  there  as  1/f  noise  with  a  spectral  density  of 
2a/7tf  arising  from  each  electron  independently,  where  a  =1/137  is  the 
fine  structure  constant.  We  will  briefly  derive  this  result  again  in  the 
present  Section,  but  we  will  stress  the  connection  between  the  coherent 
quantum  1/f  noise  and  the  usual  quantum  1/f  effect. 

Coherent  Quantum  1/f  Noise 

The  coherent  quantum  1/f  noise  will  be  derived  again  in  three  steps: 
first  we  consider  just  a  single  mode  of  the  electromagnetic  field  in  a 
coherent  state  and  calculate  the  autocorrelation  function  of  the 
fluctuations  which  arise  from  its  nonstationarity.  Then  we  calculate  the 
amplitude  with  which  this  mode  is  represented  in  the  field  of  an  electron. 
Finally,  we  take  the  product  of  the  autocorrelation  functions  calculated 
for  all  modes  with  the  amplitudes  found  in  the  previous  step. 

Let  a  mode  of  the  electromagnetic  field  be  characterized  by  the 
wave  vector  q,  the  angular  frequency  co  =  cq  and  the  polarization  X. 
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Denoting  the  variables  q  and  X  simply  by  q  in  the  labels  of  the  states,  we 
write  the  coherent  state28*29*11  of  amplitude  |Zq|  and  phase  arg  Zq  in  the 

form 


|zq>  -  exp[-(1/2)(Zq|2]  exp[Zqaq+]  |0> 


-  exp[-(1/2)|zq|2]  £  (zqn)/n!  |n>.  (68) 

n*0 

Let  us  use  a  representation  of  the  energy  eigenstates  in  terms  of  Hermite 
polynomials  Hn(x) 

|n>  -  (2nn!  Vn)-1/2  exp[-x2/2]  Hn(x)  e'ncot.  (69) 

This  yields  for  the  coherent  state  |zq>  the  representation 


Vq(x)  -  exp[-(1/2)|zq|2]exp[-x2/2]X([Zqei“‘)n/[n!(2n1((o)]1/2)Hn(x) 

n-0 

-  exp[-(1/2)|zq|2]exp[-x2/2]exp[-z2e-2io)t  +  2xze'®t].  (70) 

In  the  last  form  the  generating  function  of  the  Hermite  polynomials  was 
used11.  The  corresponding  autocorrelation  function  of  the  probability 
density  function,  obtained  by  averaging  over  the  time  t  or  the  phase  of  Zq, 
is,  for  |Zq|«1 , 

Pq(t.x)  -  <|Vqll2  I  Vq|,„2  > 

<■  (1  +  8x2|zq|2[i  +  cos  on]  -  2|Zq|2)exp[-x2/2],  (71) 

Integrating  over  x  from  ■«  to«  ,  we  find  the  autocorrelation  function 

A1  (x)  -  (2)-1/2(-j  +  2|zq|2cos  cox}  (72) 

This  result  shows  that  the  probability  contains  a  constant  background 
with  small  superposed  oscillations  of  frequency  co.  Physically,  the  small 
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oscillations  in  the  total  probability  describe  a  particle  which  has  been 
emitted,  or  created,  with  a  slightly  oscillating  rate,  and  which  is  more 
likely  to  be  found  in  a  measurement  at  a  certain  time  than  at  other  times 
in  the  same  place.  Note  that  for  Zq  -  0  the  coherent  state  becomes  the 

ground  state  of  the  oscillator  which  is  also  an  energy  eigenstate,  and 
therefore  stationary  and  free  of  oscillations. 

We  now  determine  the  amplitude  Zq  with  which  the  field  mode  q  is 

represented  in  the  physical  electron.  One  way  to  do  this11  is  to  let  a  bare 
particle  dress  itself  through  its  interaction  with  the  electromagnetic 
field,  i.e.  by  performing  first  order  perturbation  theory  with  the 
interaction  Hamiltonian 

H'  «  Aji  j4  «  -(e/c)v*A  +  ,  (73) 

where  A  is  the  vector  potential  and  <J>  the  scalar  electric  potential. 

Another  way  is  to  Fourier  expand  the  electric  potential  e/47cr  of  a  charged 
particle  in  a  box  of  volume  V.  In  both  ways  we  obtain1 1 

Izql2  >  (e/q)2(hcqV)-1 .  (74) 

Considering  now  all  modes  of  the  electromagnetic  field,  we  obtain  from 
the  single  -  mode  result  of  Eq.  (72) 

A(x)  -  C  nq{1  +  2|Zq|2cos  c»qx}  -  C{1  +  2|zq|2cos  coqx} 

-  C{1  +  2(V/23jc3)  Jd3q  |Zq|2cos  a>qx}  (75) 

Here  we  have  again  used  the  smallness  of  Zq  and  we  have  intoduced  a 
constant  C.  Using  Eq.  (74)  we  obtain 


A(  x  )  -  C{1  +  2(V/237t3)(47c/V)(e2/2hc)  J(dq/q)cos  coqx} 

■  C{1  +  2(a/ji)  Jcos(cox)dco/a>.  (76) 
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Here  a  »  e2/47ihc  is  the  fine  structure  constant  1/137.  The  first  term  in 
curly  brackets  is  unity  and  represents  the  constant  background,  or  the  d.c. 
part.  The  autocorrelation  function  for  the  relative,  or  fractional  density 
fluctuations,  or  for  current  density  fluctuations  in  the  beam  of  charged 
particles  is  obtained  therefore  by  dividing  the  second  term  in  curly 
brackets  by  the  first  term.  The  constant  C  drops  out  when  the  fractional 
fluctuations  are  considered.  According  to  the  Wiener-Khintchine  theorem, 
the  coefficient  of  cos  is  the  spectral  density  of  the  fluctuations,  S|y|2,  or 

Sj  for  the  current  density  j=e(k/m)|y|2 

S|¥|2<|v|-2>  ,  Sj<j>-2  -  2{a/jifN)  -  4.6-10-3  HN-1  (77) 

Here  we  have  included  the  total  number  N  of  charged  particles  which  are 
observed  simultaneously  in  th«  denominator,  because  the  noise 
contributions  from  each  particle  are  independent.  This  result  is  related  to 
the  well  known  conventional  Quantum  1/f  Effect1-2.  If  a  beam  of  charged 
particles  is  scattered,  passes  from  one  medium  into  another  medium  (e.g. 
at  contacts),  is  emitted,  or  is  involved  in  any  kind  of  transitions,  the 
amplitudes  zq  which  describe  its  field  will  change.  Then,  even  if  the 

initial  state  was  prepared  to  have  a  well-determined  energy,  the  final 
state  will  have  an  indefinite  energy,  with  an  uncertainty  determined  by 
the  difference  between  the  new  and  old  zq  amplitudes.  Azq.  This, 
however,  is  just  the  bremsstrahlung  amplitude  Azq.  We  thus  regain  the 
familiar  quantum  1/f  effect,  according  to  which  the  small  energy  losses 
from  bremsstahlung  of  infraquanta  yield  a  final  state  of  indefinite  energy, 
and  therefore  lead  to  fluctuations  of  the  process  rate,  or  cross  section,  of 
the  process  in  which  the  electrons  have  participated,  and  which  has 
occasioned  the  bremsstrahlung  in  the  first  place.  The  calculation  of 
piezoelectric  1/f  noise30  which  deals  with  phonons  as  infraquanta,  was 
phrased  in  terms  of  the  coherent  field  amplitudes  zq  for  the  first  time, 

although  it  is  concerned  only  with  the  usual  quantum  1/f  effect.  It  has  a 
substituted  by  the  piezoelectric  coupling  constant  g. 
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Connection  with  the  Usual  Quantum  1/f  Effect 

The  assumptions  included  in  the  derivation  of  the  above  coherent 
quantum  1/f  noise  result  are  : 

1  -  The  "bare  particle"  does  not  have  compensating  energy  fluctuations 
which  could  cancel  the  fluctuations  present  in  the  field.  The  latter  are 
due  to  the  interaction  with  distant  charges,  and  have  nothing  to  do  with 
the  bare  particle.  Therefore,  this  assumption  is  quite  reasonable. 

2  -  The  experimental  conditions  do  not  alter  the  physical  definition  of 
the  charged  particle  as  a  bare  particle  dressed  by  a  coherent  state  field. 
This  second  assumption  depends  on  the  experimental  conditions. 

One  way  to  understand  this  second  assumption  is  based  on  the 
spatial  extent  of  the  beam  of  particles  or  of  the  physical  sample 
containing  charged  particles,  and  is  specifically  based  on  the  number  of 
particles  per  unit  length  of  the  sample.  According  to  this  model,  the 
coherent  state  in  a  conductor  or  semiconductor  sample  is  the  result  of  the 
experimental  efforts  directed  towards  establishing  a  steady  and  constant 
current,  and  is  therefore  the  state  defined  by  the  collective  motion,  i.e.  by 
the  drift  of  the  current  carriers.  It  is  expressed  in  the  Hamiltonian  by  the 
magnetic  energy  Em,  per  unit  length,  of  the  current  carried  by  the  sample. 

In  very  small  samples  or  electronic  devices,  this  magnetic  energy 

Em  -  J(B2/87t)d3x  -  [nevS/c]2|n(R/r)  (78) 

is  much  smaller  than  the  total  kinetic  energy  E|<  of  the  drift  motion  of  the 
individual  carriers 

E|<  -  mv2/2  -  nSmv2/2  -  Em/s.  (79) 

Here  we  have  introduced  the  magnetic  field  B,  the  carrier  concentration  n, 
the  cross  sectional  area  S  and  radius  r  of  the  sample,  the  radius  R  of  the 
electric  circuit,  and  the  "coherent  ratio" 


s  -  Em/k  ■  2ne2s/mc2|n  (R/r)  »  2e2N7mc2, 


(80) 
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where  N'  -  nS  is  the  number  of  carriers  per  unit  length  of  the  sample  and 
the  natural  logarithm  ln(R/r)  has  been  approximated  by  one  in  the  last 
form.  We  expect  the  observed  spectral  density  of  the  mobility 
fluctuations  to  be  given  by  a  relation  of  the  form 

(1/p2)Sji(f)  -  [1/(1+s)][2ccA/fN]  +  [s/(1+s)][2cc/*fN]  (81) 

which  can  be  interpreted  as  an  expression  of  the  effective  Hooge  constant 
if  the  number  N  of  carriers  in  the  (homogeneous)  sample  is  brought  to  the 
numerator  of  the  left  hand  side.  Eq.  (81)  needs  to  be  tested 
experimentally.  In  this  equation  aA-2a(Av/c)2/37i  is  the  usual 
nonrelativistic  expression  of  the  infrared  exponent,  present  in  the 
familiar  form  of  the  quantum  1/f  effect  [4-9].  This  equation  does  not 
include  the  quantum  1/f  noise  in  the  surface  and  bulk  recombination  cross 
sections,  in  the  surface  and  bulk  trapping  centers,  in  tunneling  and 
injection  processes,  in  emission  or  in  transitions  between  two  solids. 

Note  that  the  coherence  ratio  s  introduced  here  equals  the  unity  for 
the  critical  value  N'  -  N"  >  2‘1012/cm.,  e.g.  for  a  cross  section  S  »  2*1 0*4 
cm2  of  the  sample  when  n  -  1016.  For  small  samples  with  N'«N"  only  the 
first  term  survives,  and  for  N'»N"  only  the  second  term  remains  in  Eq. 

(81).  We  hope  that  an  expression  similar  to  Eq.  (81)  will  allow  us  to 
extend  the  present  good  agreement  between  theory  and  experiment  to  the 
semiconductor  samples  of  intermediate  size. 

11.10  COHERENT  STATE  PIEZOELECTRIC  QUANTUM  1/f  NOISE 

Coherent  state  quantum  1/f  noise  has  become  known11-12  as  the 
fluctuation  caused  by  the  nonstationary  character  of  electrically  charged 
particles  or  current  carriers  in  solids.  This  nonstationarity  appears 
because  the  coherent  state  of  any  harmonic  oscillator,  in  particular  also 
the  state  of  a  physical  charged  particle  whose  electromanetic  field  is 
known  to  be  in  a  coherent  state,  is  not  an  eigenstate  of  the  Hamiltonian, 
but  rather  an  eigenstate  of  the  field  operator,  with  well  defined 
amplitude,  but  with  an  uncertain  photon  number  (see  Sec.  11.9). 
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On  the  other  hand,  the  similarity  of  the  electron  -  photon  interaction 
to  the  piezoelectric  interaction  between  electrons  and  transversal 
photons  has  led  us30  to  apply  the  conventional  quantum  1/f  theory  to  the 
case  of  piezoelectric  coupling.  The  result  of  that  calculation  is  a  rather 
large  quantum  1/f  noise  which  we  call  piezoelectric  quantum  1/f  noise,  in 
qualitative  agreement  with  the  very  large  values  given  by  the 
measurements31  >32. 

It  is  therefore  worthwhile  to  investigate  the  applicability  of  the 
concept  of  coherent  states  1/f  noise  to  the  case  of  transversal  phonons  as 
infraquanta.  In  this  Section  we  calculate  for  the  first  time  the  coherent 
state  quantum  1/f  noise  caused  by  phonons  as  infraquanta. 

Starting  from  the  expression30 

Aq(k)  =*  (2;ig/qV)1/2  (vs/q)  [hq/2m  -  hq0-k/m  +vs]~ 1  (82) 

of  the  amplitude  of  the  phonon  state  q  associated  with  the  coherent  state 
which  describes  the  deformation  field  of  an  electron  of  wave  vector  k,  we 
obtain  the  number  of  phonons  emitted  when  the  state  was  created,  by 
squaring  this  amplitude  and  summing  over  the  phonon  modes: 

Nphcoh  -  (2)[V/(2n)3]  Jd3q  (2jtgvs2/Vq3)[vs  -  hq„.k/m]-2 

-  (2)  (g/it)  [1-(hk/mvs)2]-l  Jdq/q.  (83) 

This  is  also  the  number  of  phonons  obtained  by  expanding  the  lattice 
strain  (or  deformation)  field  of  the  electron  in  terms  of  phonon  modes. 
Here  g  is  the  piezoelectric  coupling  constant,  V  the  volume  of  a 
normalization  box  which  disappears  later  in  the  calculation,  vs  the  speed 
of  sound,  q0  *  q/q  a  unit  vector,  and  the  factor  2  placed  in  brackets 
accounts  for  the  two  polarization  modes  associated  with  each  phonon 
vector  q.  The  angular  integration  was  performed  with  the  assumption 
hk/mvs  <  1,  i.e.,  for  subsonic  electrons.  For  higher  momenta  no  1/q 
spectrum  is  obtained,  and  so  no  1/f  contributions  are  added.  This  implies 
a  decrease  of  the  1/f  noise  at  supersonic  drift  velocities  of  the  carriers, 
compared  to  the  extrapolation  from  the  subsonic  region. 
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According  to  the  general  quantum  1/f  theoretical  relation2,  the 
phonon  number  spectral  density  given  by  the  last  equation  is  also  half  the 
spectral  density  of  fractional  coherent  state  quantum  1/f  current 
fluctuations  observable  in  a  beam  of  electrons  of  wave  vector  k: 

Sj(f)/j2  -  4  (g/f)  [1-(bk/mvs)2]-i.  (84) 

This  expression  differs  from  our  previous  (piezoelectric  conventional 
quantum  1/f)  result  by  the  absence  of  the  factor  (1/3)[h(k  -  k')/mvs]2. 

A  more  detailed  presentation  of  this  derivation  is  found  in  a  recently 
published  article33 


III.  APPLICATIONS  OF  THE  QUANTUM  1/F  THEORY 

111.1  DERIVATION  OF  MOBILITY  QUANTUM  1/f  NOISE  IN  n+-p  DIODES 

For  a  diffusion  limited  n+-p  junction  the  current  is  controlled  by 
diffusion  of  electrons  into  the  p  -  region  over  a  distance  of  the  order  of 
the  diffusion  length  L  -  (Dptn)1^  which  is  shorter  than  the  length  wp  of 

the  p  -  region  in  the  case  of  a  long  diode.  If  N(x)  is  the  number  of 
electrons  per  unit  length  and  Dn  their  diffusion  constant,  the  electron 

current  at  x  is 


Ind  *  "  ©DpdN/dx,  (85) 

where  we  have  assumed  a  planar  junction  and  taken  the  origin  x  *  0  in  the 
junction  plane.  Diffusion  constant  fluctuations,  given  by  kT/e  times  the 
mobility  fluctuations,  will  lead  to  local  current  fluctuations  in  the 
interval  Ax 

8Alnc|(x,t)  *  lndAxSDn(x,t)/Dn.  (86) 

The  normalized  weight  with  which  these  local  fluctuations  representative 
of  the  interval  Ax  contribute  to  the  total  current  lq  through  the  diode  at  x 


43 


■  0  is  determined  by  the  appropriate  Green  function  and  can  be  shown  to  be 
(1/L)exp(-x/L)  for  Wp/L  »  1.  Therefore  the  contribution  of  the  section  Ax 

is 

5Al<j(x,t)  -  (Ax/L)exp(-x/L)Jncj5Dn(x,t)/Dn,  (87) 

with  the  spectral  density 

SAId(x.f)  -  (AX/L>2  exp(-2x/L)  ln<j2  SD„(x,f)/Dn2.  (88) 

For  mobility  and  diffusion  fluctuations  the  fractional  spectral  density  is 
given  by  annd^NAx,  where  annd  is  determined  from  quantum  1/f  theory 
according  to  Sec  111.3.  With  Eq.  (85)  we  obtain  then 

SA|d(*.f)  -  (Ax/L2)  exp(-2x/L)  (eDndN/dx)2  <XHnd/fN-  (89) 

The  electrons  are  distributed  according  to  the  solution  of  the  diffusion 
equation,  i.e. 

N(x)  »  [N(0)  -  Np]  exp(-x/L);  dN/dx  *  -{[N(0)  -  Np]/L}  exp(-x/L).  (90) 

Substituting  into  Eq.  (89)  and  simply  summing  over  the  uncorrelated 
contributions  of  all  intervals  Ax  we  obtain 

wp 

S|d(0  -  oHnd(eDn/L2)2  J[N(0)  -  Np]2e-4x/Ldx  /{[N(0)  .  Np]  e-x/L  +  Np}.  (9i) 

0 

We  note  that  eDn/L2  -  e/xn-  With  the  expression  of  the  saturation  current 
Ip  -  e(Dn/xn)1/2Np  and  of  the  current  I  -  lo[exp(eV/kT)  -  1],  we  can  carry 
out  the  integration 


S|d (0  -  aHnd(el^fxn)  Ja2u3du/(au  +  1)  -  aHnd(e,/fxn)  F(a)-  (92) 


Here  we  have  introduced  the  notations 
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u  -  exp(-x/L),  a  -  exp(eV/kT)  -  1 , 

F(a)  -  1/3  -  1/2a  +  1/a2  -  (1/a3)ln(1+a).  (93) 

Eq.  (92),  obtained  by  van  der  Ziel  and  Anderson,  gives  the  diffusion  noise 
as  a  function  of  the  quantum  1/f  noise  parameter  aund-  A  similar  result 
can  be  derived  for  the  quantum  1/f  fluctuations  of  the  recombination  rate 
in  the  bulk  of  the  p  -  region. 

111.2.  OVERVIEW  OF  QUANTUM  1/f  NOISE  IN  n+p  JUNCTIONS  AND  MIS  DIODES 

As  we  have  seen  in  the  Sec.  II,  quantum  1/f  noise  is  a  low-frequency 
fluctuation  process  present  in  elementary  cross  sections  o  and  process 
rates  r,  given  by  the  fractional  spectral  density 

Sa(f)/a2  -  Sr(f)/T2  -  (4a/3jcfN)(Av/c)2  (94) 

for  conventional  quantum  1/f  noise  which  is  applicable  to  small  devices, 
and 


Sa(f)/o2  -  Sr(f)/T2  -  (2a/7cfN)  -  ac/fN  (95) 

for  coherent  state  quantum  1/f  noise  which  is  applicable  to  large  devices 
in  which  the  energy  of  the  carrier  drift  motion  is  predominantly  magnetic, 
rather  than  kinetic.  Here  Av  is  the  velocity  change  of  the  carriers  in  the 
processes  considered,  a  -  1/137  is  the  fine  structure  constant,  N  the 
number  of  carriers  simultaneously  interrogating  the  cross  sections  or 
process  rates,  ac  -  4.6  10'3  the  coherent  quantum  1/f  noise  coefficient, 
and  c  the  speed  of  light.  The  two  forms  of  quantum  1/f  noise  are  closely 
related  infrared  divergence  phenomena  which  arise  due  to  the  interaction 
of  electrons  and  soft  photons,  as  we  have  seen  in  Sec.  11.9. 

Both  in  n+p  diodes13*16'18-23-3*  anc|  metal-insulator-semiconduc- 
tor35-38  (MIS)  devices  the  current  will  be  determined  by  cross  sections  crj 
such  as  recombination  and  scattering  cross  sections  (by  phonons  and 
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lattice  defects),  as  well  as  by  other  process  rates  (e.g.,  band  to  band  and 
trap-assisted  tunneling,  particularly  important  in  long-wavelength  Hgi_x 
CdxTe  detectors).  The  spectral  density  of  the  resulting  quantum  1/f 
current  fluctuations  is  therefore 

S,(f)  -  Si  [ai/3oi]2Sai(f)  +  S  [3l/ari]2Sri(f),  (96) 

where  the  spectral  densities  on  the  right  hand  side  will  be  given  by  Eqs. 
(94)  or  (95),  depending  on  the  size  of  the  device. 

A.  N+P  Junctions 

The  current  density  I  of  Eq.  (96)  contains  a  diffusion  term  ld,  a  term 
lr  caused  by  recombination  in  the  space  charge  region,  a  surface 
recombination  term  ls,  a  tunneling  term  lt  and  a  photovoltaic  term  caused 
by  the  creation  of  electron  hole  pairs  by  photons: 

I  »  Id  +  lr  +  Is  +  *t  +  OH  ^ 

-  qni{(ni/n0)(nn/Tn)i/2(eqV/kT.1)  +(W/t)(eqV/2kT.-|)  +S}  +|t  +qri<J>.  (97) 

Here  n;  is  the  intrinsic  concentration,  n0  the  concentration  of  acceptors  on 
the  p  side,  Dn  and  xn  the  diffusion  constant  and  lifetime  of  minority 
carriers  on  the  p  side,  W  the  width  of  the  depletion  region,  x  *  xpo  +  no  the 
Shockley-Hall-Read  lifetime,  V  the  applied  voltage,  s  the  surface 
recombination  speed,  i\  the  quantum  efficiency  and  4>  the  incident  flux  of 
photons.  With  the  exception  of  the  last  term,  the  terms  in  Eq.  (97)  are 
known  as  dark  current  components. 

The  first  term  in  rectangular  brackets  in  Eq.  (97)  gives  the  diffusion 
current  density  Id,  and  yields  a  noise  term16-18 

S,d(f)  -  (ad+ar)(qld/fxn)F(a);  ad  ■  (4a/37i)(h/m*bc)2exp(-e/2T)  (98) 

in  Eq.  (96)  for  small  devices,  and 
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S|d(f)  -  ac(qld/fxn)F(a)  (99) 

for  large  devices,  with  F(a)  -  -1/3  -1/2a  +1/a2  -(1/a3)ln(1+a)  and  a  - 
exp(qV/kT)  -1.  Here  b  is  the  lattice  constant,  m*  the  effective  mass  of 
the  electrons,  and  6  the  Debye  temperature.  These  results  are  obtained  by 
adding  the  scattering  and  recombination  quantum  1/f  cross  section 
fluctuation  spectra  for  all  points  on  the  p  side  of  a  long  n+p  diode  with  the 
appropriate  Green  function  weight  [3l/3aj]2,  as  prescribed  by  Eq.  (96).  The 
spatial  dependence  exp(-x/Ld)  of  this  Green  function  is  given  by  the 
diffusion  equation  applied  to  electrons  on  the  p  side,  with  being  the 
electron  diffusion  length.  If  the  length  wp  of  the  p  side  is  shorter  than  !_<*, 
we  have  to  replace  F(a)  by  F(a)  -  F[a  exp(-wp/Ld)]  in  Eqs.  (98)  -  (99),  and 
also  to  include  the  quantum  1/f  noise  of  the  recombination  cross  sections 
at  the  p  contact. 

Similarly  we  obtain16*20 

S|r  -  arqlrtanh(qV/2kT)/fx; 

-  (4cc/37cc2)[2q(Vdiff  -V)  +3kTJ{(mn*)1/2  +(mp*)1/2}-2.  (100) 

The  corresponding  expression  for  the  surface  recombination  current  is 
given  by  Eq.  (100)  with  half  of  the  surface  potential  jump  U  added  to  the 
diffusion  potential  Vdiff.  The  tunneling  current  It  and  its  noise  is  the  same 
as  in  MIS  devices, and  will  be  discussed  below.  In  the  final  form  of  Eq. 

(96)  the  fine  structure  constant  a  appears  as  a  general  factor. 

2.  MIS  Devices 

The  MIS  infrared  detector  devices35  considered  here  are  homogeneous 
Hgt.xCdxTe  slabs  of  thickness  of  the  order  of  the  minority  carrier 
diffusion  length,  on  which  insulated  metallic  field  effect  plates  are 
applied  with  a  potential  that  creates  a  surface  inversion  layer.  The 
transition  between  the  bulk  and  the  inverted  surface  layer  is  a  pn  junction 
which  is  similar  to  a  conventional  pn  junction  caused  by  inhomogeneous 
doping,  connected  in  series  with  a  capacitor,  but  is  free  of  the  defects 
introduced  in  the  doping  process.  The  theory  summarized  for  junctions  in 
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the  previous  Section  is  again  applicable  for  the  quantum  1/f  noise  of  the 
(dark)  current  flowing  through  the  junction,  but  this  time  the  longer  life 
times  in  the  denominator  of  Eqs.  (98),  (99)  and  (100)  will  lead  to  lower 
quantum  1/f  noise.  This  is  a  considerable  advantage  of  MIS  devices,  but  on 
the  other  hand,  they  have  to  be  operated  in  a  pulsed  mode  in  which  the 
photoelectrically  generated  carriers  separated  by  the  field  of  the  induced 
junction,  and  accumulated  in  the  potential  well  of  the  surface  inversion 
layer,  have  to  be  removed  and  registered  periodically,  by  alternating  the 
inverting  voltage  pulses  at  the  plate  with  short  "readout"  potentials  which 
generate  temporary  conditions  close  to  flat  energy  bands  at  the  surface. 

In  addition  to  the  photoelectrically  generated  carriers,  the  readout 
signal  also  contains  the  carriers  accumulated  through  the  various  dark 
current  mechanisms  mentioned  above.  Therefore,  at  low  frequencies  the 
1/f  noise  contained  in  the  dark  current  I  will  limit  the  performance  of  MIS 
detectors  as  infrared  detectors. 

The  integration  implied  by  the  pulsed  mode  of  operation  acts  as  a  low 
pass  filter  which  does  not  influence  the  signal  to  noise  ratio  at 
frequencies  well  below  the  readout  frequency.  Furthermore,  the  pulsed 
mode  of  operation  does  not  affect  the  quantum  1/f  noise,  because  the 
quantum  1/f  fluctuations  of  cross  sections  and  process  rates  are 
continuously  tested  by  carriers  in  thermal  agitation,  and  are  present 
independently  of  the  method  of  sampling  used  by  us. 

In  general  the  dark  current  I  can  again  be  written  in  the  form 


I  =*  ld  +  lr  +  ls  +  It  +  (101) 

The  five  terms  on  the  right  hand  side  correspond  to  minority  carrier 
diffusion  from  the  bulk,  generation  from  Shockley-Read  centers  in  the 
depletion  region,  generation  from  SR  centers  at  the  surface,  tunneling,  and 
photoelectric  generation  by  the  thermal  radiation  background  flux  o.  We 
shall  give  the  formulae20  derived  in  detail  in  Sec.  111.4,  which  determine 
these  terms  below,  also  including  an  example  of  their  calculation  in  a  p- 
type  device 
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Id  -  (qnj2/p0)[kWqtn]1/2  -  (1.6  lO-^C  36  1024  cm-6/10i5Cm3)  [O.OIV 


1.5  105(cm2/Vs)/10-6s]1/2  «  2  10-4A/cm2,  (102) 

lr  -  qniW/2-o  *  [(1.6  lO-^C  6  1012cm-3  2  10-4cm  /(2  10-6s)] 

»  10*4A/cm2  (103) 

ls  -  Jb  -  qrijS/2  -  [(1.6  10-19C  6  lO^cnr3  20cm/s] 

-  1.8  10'5A/cm2,  (104) 

lt*10-l3(NrOs/Eg)exp[-5.3  10<5Eg2/E]-1 0-13(4  10l2Cm-30.2V/0.062V) 

exp[-5.3  106(0.062V)2/3000V/cm]  »  8  lO^A/cmS,  (105) 


Here  n;  -  6  1012cm3  is  the  intrinsic  carrier  concentration,  and  the 
concentration  of  holes  was  taken  to  be  p0  -  1015cnr3.  The  mobility  p  = 

1.5  1015cm2/Vs  as  well  as  the  life  time  xn  ■  10-6s  of  the  minority 
carriers  have  to  be  replaced  by  pp  and  xp  for  the  case  of  n-type  devices. 

The  surface  recombination  speed  was  denoted  by  s  *  40cm/s,  and  the 
concentration  of  intermediary  states  effective  for  tunneling  was  denoted 
by  Nr  -  4  10i2cm-3.  The  band  gap  considered  was  Eg  *  0.065eV,  the 

surface  potential  <J>S  -  0.2V  and  the  electric  field  below  the  surface  E  = 
3000V/cm.  All  numerical  values  have  been  included  only  as  an  example 
and  are  not  characteristic  of  a  particular  device.  The  numerical  factors  in 
Eq.  (105)  correspond  to  p  -  type  HgCdTe  with  a  20  pm  cutoff  wavelength 
and  were  taken  from  the  paper  of  Kinch  and  Beck36.  For  n-type  devices  we 
also  need  to  include  tunneling  via  surface  states.  If  the  density  of  fast 
surface  states  is  denoted  by  NfS  *  1012  cm-zv*1,  the  current  generated 
from  tunneling  via  a  uniform  density  of  fast  surface  states  across  the 
band  gap  will  be  given  by3  5 

•tsc  -  -qPiNfskT/2q  -  1.75  10*4  A/cm2.  (106) 

The  sum  of  the  first  four  currents  on  the  right  hand  side  of  Eq.  (101) 
must  be  smaller  than  the  fifth  which  corresponds  to  the  thermal  radiation 
background  flux,  for  background  limited  (BLIP)  performance.  Although  not 
all  terms  in  the  dark  current  are  of  importance,  we  still  retain  them, 
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because  their  quantum  1/f  noise  may  be  quite  significant,  even  if  the 
corresponding  current  is  negligible.  We  shall  now  proceed  with  the 
calculation  of  quantum  1/f  noise  contributions  from  all  these  currents. 

The  quantum  1/f  noise  in  Id  and  lr  is  the  same  as  in  n+p  diodes,  and  is 
given  by  Eqs.  (98)-(1 00). 

The  surface  recombination  current  is  a  dark  current  component 
originating  from  surface  states  at  the  interface  between  the  surface 
passivation  layer  and  the  bulk.  The  recombination  cross  sections  of  these 
states  will  exhibit  quantum  1/f  noise.  The  quantum  1/f  noise  coefficient 
an  reflects  the  velocity  change  of  carriers  involved  in  this  process 

<xH  -  (4ct/37t)(2/m*c2)[3kT/2  +  ell/2  +  0.1  eV],  (107) 

where  U  is  the  surface  potential  jump  present  between  the  surface 
passivation  layer  and  the  bulk  even  if  no  gate  voltage  is  applied,  and  V  is 
the  applied  gate  voltage  which  we  take  with  a  coefficient  less  than  unity, 
here  for  example  with  a  weight  of  0.1. 

The  calculation  of  surface  recombination  quantum  1/f  noise  is 
similar  to  the  calculation  of  quantum  1/f  noise  from  recombination  in  the 
space  charge  region.  However,  in  this  case  the  cross  sections  are  not 
distributed  over  the  width  of  the  junction,  but  rather  are  concentrated  at 
the  surface  which  is  characterized  by  the  surface  potential  <X>S.  Therefore, 
with  x=eV/2kT,  we  can  write  an  expression  of  the  form 

S|S(f)  »  aHelstanh  x  /f(xno  +  xpo)exp(eOs/kT).  (108) 

We  note  that  the  fabrication  process  of  MIS  structures  introduces 
less  bulk  defects  than  th6  fabrication  of  photovoltaic  devices. 

Nevertheless,  the  fabrication  of  MIS  devices  introduces  some  defects  in 
the  bulk  layer  located  right  under  the  surface.  These  defects  will 
manifest  themselves  through  a  contribution  to  indirect  tunneling. 

Quantum  1/f  Noise  in  the  Tunneling  Rate.  The  tunneling  component  It 
of  the  dark  current  is  particularly  important  in  narrow  band  gap  material. 
In  the  case  of  tunneling  from  a  surface  accumulation  layer  to  the  bulk,  the 
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velocity  change  of  the  carriers  will  lead  us  from  the  thermal  velocity  of 
the  carrier  on  one  side  of  the  barrier  to  the  thermal  velocity  of  a  carrier 
on  the  other  side  of  the  barrier,  if  band  to  band  tunneling  is  considered. 

If,  however,  the  tunneling  goes  via  intermediate  states  located  in  the  band 
gap,  the  velocity  is  zero  as  long  as  the  carrier  is  stationary  in  the 
intermediate  state.  We  can  therefore  write  the  1/f  noise  coefficient 

<xH  -  (4a/37t)6kT/m*c2  (109) 

for  band  to  band  tunneling,  and 

an  -  (4a/37t)3kT/m*c2  (110) 

for  tunneling  via  intermediate  states  in  the  band  gap,  where  we  have 
considered  the  average  squared  velocity  change  two  times  smaller.  The 
effective  mass  is  the  mass  of  the  minority  carriers  in  the  bulk  material. 

For  band  to  band  tunneling  from  a  surface  inversion  layer  to  the  bulk, 
the  velocity  change  of  carriers  corresponds  to  an  energy  difference  of  the 
order  of  the  band  gap  Eg  plus  an  energy  difference  of  the  order  of  the 
thermal  energy  3kT/2,  provided  we  are  dealing  with  deep  inversion,  as 
used  in  practical  MIS  devices  in  the  pulsed  mode  of  operation.  This  yields 
the  quantum  1/f  coefficient 

aH  »  (4a/3rc)(Eg  +  3kT/2)/m*c2  (111) 

for  band  to  band  tunneling.  For  tunneling  via  intermediate  states  in  the 

band  gap  the  corresponding  energy  difference  will  be  smaller,  and 
therefore  we  replace  Eg  by  Eg/2 

aH  -  (4a/37c)(Eg/2  +  3kT/2)m*c2  (112) 

This  relation  is  applicable  both  if  the  intermediate  states  are  located  in 
the  depletion  region  or  at  the  surface. 

In  the  last  four  equations  we  did  not  divide  by  the  number  of  carriers 
simultaneously  involved  in  the  tunneling  process,  because  this  number  is 
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less  than  unity  for  practical  tunneling  currents.  Whenever  a  cross  section 
or  a  process  rate  is  tested  with  one  electron  or  less  than  one  electron  at  a 
time,  the  effective  number  of  electrons  in  the  denominator  of  the  quantum 
1/f  formula  must  be  replaced  by  unity.  Let  us  calculate  the  average 
number  of  carriers  simultaneously  present  in  the  tunneling  process  at  any 
time.  The  tunneling  process  occurs  over  a  distance  d  -  Eg/eE,  and  the 
speed  v  of  the  carriers  will  be  of  the  order  of  the  thermal  speed  in  the 
case  of  an  accumulation  layer,  and  of  the  order  of  the  band  gap  energy  in 
the  case  of  a  surface  inversion  layer.  Dividing  d  by  v,  we  obtain  a 
tunneling  time  t  =*  6  10'13s  for  accumulation  layers  and  t  =*  3  10'13s  for 
inversion  layers.  From  Eq.  (105)  we  know  that  the  tunneling  current  is  of 
the  order  of  10*3A/cm2.  Multiplying  this  by  t,  we  obtain  3  -  6  10-16  C, 
i.e.,  2000  -  3000  electrons/cm2  tunneling  simultaneously  in  a  device  of  1 
cm2  area.  In  a  device  of  dimensions  50  jim  x  50  pm  «  2.5  10‘5  cm2  the 
average  number  of  carriers  in  the  process  of  tunneling  at  any  time  is 
therefore  0.05  -  0.075,  and  this  is  indeed  much  less  than  unity. 
Nevertheless,  if  the  area  of  the  device  exceeds  5  10*4  cm2,  Eqs.  (109)  - 
(112)  require  an  additional  factor  e/tltA  which  makes  the  noise  spectral 
density  proportional  to  lt  and  A,  rather  than  to  the  square  of  these 
quantities. 

The  photoelectric  current  will  reflect  the  fluctuations  in  the  number 
of  photons  arriving  from  the  radiation  background.  The  quantum  efficiency 
will  not  exhibit  considerable  quantum  1/f  noise,  because  the  generated 
carriers  will  be  collected  with  certainty.  Therefore  the  collection  of 
photoelectrically  generated  carriers  is  not  controlled  by  any  cross  section 
or  process  rate  affected  by  considerable  quantum  1/f  noise.  If  we  neglect 
the  1/f  noise  generated  in  the  series  resistance  of  the  diode,  there  should 
be  no  photoelectrically  generated  1/f  noise  from  a  short-circuited  diode. 

The  resulting  1/f  noise.  Since  the  various  dark  current  fluctuations 
with  1/f  spectrum  are  statistically  independent,  the  total  1/f  noise  is 
simply  obtained  by  summing  all  contributions. 

Let  us  now  evaluate  the  spectral  density  S'(f)  of  fractional 
fluctuations  in  the  various  dark  current  noise  contributions  per  cm2  of 
transversal  detector  (or  gate)  area,  including  also  a  numerical  example  for 
a  MIS  infrared  detector.  For  a  given  detector,  this  needs  to  be  divided  by 
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the  area  A  of  the  detector  to  yield  the  corresponding  fractional  spectral 
densities:  S(f)  =[S’(f)]/A.  Fractional  fluctuations  are  dimensionless,  so 
S'(f)  will  have  the  dimension  of  a  reciprocal  frequency,  times  a  squared 
length  unit  which  simplifies  when  we  divide  by  the  area  of  the  detector  at 
hand.  Let  S'|d  be  the  spectral  density  of  fractional  fluctuations  in  the 
noise  caused  by  quantum  1/f  fluctuations  in  diffusion,  S'\T  in  bulk 
recombination,  S'|s  in  surface  recombination,  and  S’|t  in  tunneling.  With 
mp*  =  0.55  m0,  mn*  *  0.02  m0,  xn  -  10'6s,  Eg  *  0.1  eV,  3kT/2  =  0.01  eV,  we 
obtain  for  a  p-type  MIS  device  with  wp>>Ld 

S’ld  *  («Hnd  +  «Hnr)[e/fxnld]F(a)  =  aco^el^/^kT^TnJ^^pJJFOVa 
=(4.6  10'3/4fNp)  4  10*10C1/2/[10'6S  1.5  105(cm2/Vs)  4  10'21  J]1/2 
=1.8  10-6cm2/f,  [or  10'8  cm2/f,  with  ap  =2  10*5  for  incoher.  noise]  (113) 

S'ir  »  [aHee/f(xno  +  xpo)lr]tanhx  =  [aH0e/(feAwn  tanhx)]tanhx 

=  aH0/fAwnj  =  4.6  10-9cm2/f  (x  =  eV/2kT)  (1 14) 

S’,s  =  (4a/37i)(2/nrc2)[3kT/2  +  eU/2  +  O.IVe]  [e  tanhx  /f(xno+xpo)ls] 

=(4a/37i0. 02) (2/500, 000) [0.025+0. 5+0. 5]  [e  tanhx/feAwn^ex-l )] 

=  7  10-8cm2/f  (115) 

S'|t  =  (4a/3it)(Eg  +  3kT/2)/m*c2  =  (4/9.5  137  0.02)(0.1 1/500,000) 

»  3.3  10-8cm2/f  (116) 

S'|d  was  calculated  in  the  small  bias  limit  for  wp»L,  but  wp=  0.25  L  gives 
the  same  result;  the  incoherent  case  with  a  lattice  constant  of  b  =  0.65 
nm  and  0  =  320  K  was  also  listed  above,  and  would  give  1.8  10'10cm2/f 
for  a  n-type  device.  Eqs.  (113)  -  (116)  would  be  reduced  mp*  /mn*  =  27.5 
times  for  n-type  devices.  S*|s  has  been  calculated  with  the  inclusion  of  a 
term  of  10%  of  the  applied  gate  voltage  V  (■  5  volt  in  our  example)  into 
the  kinetic  energy  of  the  carriers  at  the  surface. 
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We  conclude  that  the  effective  mass  present  in  the  denominator 
leads  to  lower  conventional  quantum  1/f  noise  in  n  type  devices. 

111.3  QUANTUM  1/F  NOISE  IN  JUNCTION  INFRARED  DETECTORS 

A.  QUANTUM  1/f  NOISE  SOURCES  APPLICABLE  TO  n+p  DIODES 

The  various  quantum  1/f  noise  sources  present  in  n+p  diodes  which  we 
will  consider  here,  have  in  common  that  they  can  be  described  by  Eq.  (3) 
with  a  Av2  value  described  by  an  effective  energy  E  =  qV0,  where  V0  is  an 
effective  potential  which  accelerates  the  carriers 

(Av2)/c2  -  2E/m*c2  -  2qV0/m*c2.  (117) 

In  various  quantum  1/f  noise  sources  the  applicable  velocity  change  v 
determines  the  Hooge  constant: 

(a)  Recombination  quantum  1/f  noise  in  the  bulk  space  charge  region 16: 

S|(f)  =»  anq|lgr|[tanh(qV/2kT)]/f(Tno  +  ^po)'  O  18) 

aH  -  (4a/3jt)[2q(Vdif-V)  +  3kT]{(mn‘)1/2  +  (mp*)1'2}-2c-2, 

where  lgr  -  qAWnj(  eQV/2kT  -  1)/(xp0  +  tpo)  is  the  recombination  current, 
xno  and  tpo  the  Shockley-Hall-Read  lifetimes,  V  the  applied  voltage  and 
V(jjf  the  diffusion  potential  of  the  junction.  Introducing  an  "effective 
carrier  number"  Neff  -  IlgrK^no+^poVqltanhfaV/2^)],  Eq.  (118)  may  be 
written  in  the  form 


S|(f)  -  aHlgr2/fNeff.  (119) 

(b)  Quantum  1/f  noise  in  the  surface  recombination  current  of  n+-p 
diodes.  This  effect  is  caused  by  quantum  1/f  fluctuations  of  the  surface 
recombination  cross  sections.  The  calculation  is  similar  to  the  previous 
(bulk)  case,  but  the  GR  process  is  localized  at  the  surface,  and  the 
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additional  electric  field  arising  from  the  potential  jump  2U  at  the 
interface  between  the  bulk  and  the  oxide  and  passivation  layers  will  lead 
to  increased  velocity  changes  of  the  carriers  in  the  recombination  process 
and  to  larger  h  values.  Including  also  the  quantum  1/f  mobility 

fluctuation  noise  of  the  spreading  resistance  caused  by  the  concentration 
of  generation  and  recombination  currents  at  the  intersection  of  the 
depletion  region  with  the  surface  of  the  diode,  we  obtain  a  current  noise 
contribution 

S|(f)  =  a'Hq|l'grl[tanh(qV/2kT)]/f(t'no  +  ^po) 

+  ap(l'gr)2 '*nfP  (ln(A/w2)]2,  (120) 

oh  =  (4a/3K)[2q(Vdif+U-V)+3kT]{(mn‘)1/2  +  (mp*)1/2}-2c-2  (121) 

Here  the  primed  quantities  refer  to  the  surface,  i.e.  we  have  introduced 
the  surface  GR  current  Pg r,  the  lifetimes  in  the  vicinity  of  the  surface  and 
the  a’n  parameter  for  surface  recombination.  P  is  the  perimeter,  A  the 

area  and  W  the  width  of  the  junction.  The  quantum  1/f  mobility 
fluctuation  part  is  expressed  in  terms  of  the  global  ctu  parameter  which 
includes  all  types  of  scattering  weighed  with  the  appropriate  mobility 
ratio  factors.  Introducing  again  an  "effective  carrier  number"  N’eff  = 
H’grK^’no+'t'poVQttanh^V/2^)],  the  first  term  of  Eq.  (120)  may  also  be 
written  in  the  form: 

S|(0  -  a'H(l'gr)2/fN'eff.  (122) 

Due  both  to  the  surface  potential  jump  2U  and  the  1/f  noise  of  the 
spreading  resistance,  the  surface  recombination  current  will  be  noisier 
than  an  equal  bulk  recombination  current,  and  this  is  in  agreement  with 
the  experimental  data34. 

(c)  Injection  -  extraction  quantum  1/f  noise16,  due  to  injection  or 
extraction  of  carriers  across  barriers.  In  this  case,  for  not  too  small 
currents 


S|(f)  -  aHniPIq/ft; 
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aHni  -  (4a/3n)[2q(Vdif-V)+3kT]/mn*c2,  (123) 

where  I  is  the  injected  current,  and  is  the  time  of  passage  of  a  carrier 
through  the  barrier  region.  Introducing  again  an  effective  carrier  number 
Neff  *  Ix/q,  Eq.  (123)  may  be  written  in  a  more  general  form,  valid  also  for 

very  small  I 


S,(f)  -  aHnil2/fNeff.  (124) 

Note  that  in  each  case  Neff,  as  long  as  it  is  larger  than  1,  is  proportional 
to  I  (otherwise  Neff  -  1,  see  below),  and  that  a h  depends  on  bias. 

(d)  Recombination  in  the  p  region  of  a  long  n+-p  diode  of  length  W  » 


(Dnxn)l/2  yields,  as  shown  in  Sec.  Ill.f, 

S|(f)  -  OHnrll|qF(a)/ftn  -  OHnrl2/fNeff;  (125) 

Neff  =  IIIWqF(a),  aHnr  -  (4a/3i)(3kT/mn*c2],  (126) 

Here  F(a)  is  given  by  Eq.  (93): 

F(a)  *  1/3  -  1/2a  +  1/a2  -  (1/a3)ln(1+a),  (127) 

where 

a  -  exp(eV/kT)  - 1 .  (128) 


In  this  case  aH  is  very  small.  Therefore  this  contribution  will  often  be 
masked  by  mobility  fluctuation  1/f  noise. 

(e)  Recombination  quantum  1/f  noise  at  the  p  -  region  contact.  The 
carriers  lose  an  average  energy  3kT/2  in  the  conduction  band  in  the 
recombination  process,  gain  an  energy  of  the  order  of  qEg  due  to  energy 
differences  of  the  order  of  the  band  gap  Eg  accelerating  the  carriers  at  the 

contact  surface,  and  gain  an  average  energy  3kT/2  in  the  valence  band,  so 
that 


«Hnc  -  (4<x/37i)[(2qEg  +  3kT)/mn*c2]. 


(1  29) 
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This  effect  will  be  important  in  short  diodes  (Wp  <  L),  where  the 
recombination  speed  at  the  contact  partially  determines  the  current 
through  the  device. 

(f)  Quantum  1/f  noise  in  the  tunneling  rate.  Tunneling  is  observed  in 
n+-p  diodes  with  sufficient  gate  bias^4.  If  we  assume  that  the  momentum 
change  of  the  carriers  in  the  tunneling  process  is  of  the  order  of  the 
thermal  r.m.s.  momentum,  we  obtain  a  minimal  quantum  1/f  noise  power 
spectrum 

S|(f)  -  aHl2/fNeff;  Neff  -  |l|t/q,  (130) 

aH  -  (4a/3jt)[3kT/mn‘c2],  (131) 

where  is  the  time  of  passage  through  the  barrier,  or  tunneling,  i.e.  the 
time  during  which  jach  carrier  contributes  to  the  current  through  the 
barrier.  Since  the  width  of  the  barrier  crossed  by  tunneling  is  small,  this 
time  is  very  short,  of  the  order  of  10‘14s.  Neff  will  then  become  larger 

than  1  at  currents  exceeding  10’5a,  leading  to  a  linear  current  dependence 
of  the  noise  power.  At  lower  bias  Neff  must  be  set  equal  to  unity  in  Eq. 

(130).  and  this  gives  a  quadratic  current  dependence.  On  the  other  hand,  at 
large  applied  voltages,  or  for  tunneling  out  of  a  metal,  kT  has  to  be 
replaced  by  a  larger  energy  in  the  expression  of  an- 

(g)  Mobility  fluctuation  1/f  noise  in  a  long  n+-p  diode.  Although 
characterized  mainly  by  a  well  defined  momentum  transfer  of  the  carriers 
rather  than  by  an  effective  energy,  and  although  we  have  discussed  it 
already  above,  this  important  process  is  listed  here  again.  As  was  shown 
in  Sec.  111.1,  this  contribution  can  be  written  in  the  form 

S|(f)  -  aH|l|qF(a)/ftn  -  oHl2/fN0ff;  Neff  -  |l|x/qF(a).  (132) 

where  an  is  given  by 


aH  -  (4a/3ji)[h/m‘ac]2exp(-  /2T)  +  (4a/37t)[6kT/mc2]  (133) 
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for  small  devices,  and  is  usually  larger  than  in  the  case  (d),  particularly  if 
the  effective  mass  of  the  carriers  is  small.  Here  is  the  Debye 
temperature,  a  the  lattice  constant,  and  T  the  temperature.  The  first 
term  in  Eq.  (133)  is  from  Umklapp  and  inter  -  valley  lattice  scattering, 
while  the  second  describes  normal  scattering  processes.  The  function  F(a) 
is  again  given  by  Eq.  (93).  In  general  this  contribution  includes  Umklapp 
scattering,  intervalley  scattering,  normal  (non  -  Umklapp)  phonon 
scattering,  impurity  scattering,  and  optical  phonon  scattering 
contributions  in  the  quantum  1/f  mobility  fluctuation  Hooge  parameter  . 
For  large  samples  or  devices,  Eq.  (95)  is  applicable.  Eq.  (81)  is  a 
suggestion  for  the  transition  region. 

B  DISCUSSION  AND  RECOMMENDATIONS 

The  quantum  1/f  noise  formulae  presented  above  have  been  applied  by 
Radford  and  Jones34  to  1/f  noise  in  GR,  diffusion  and  tunneling  currents  in 
both  double  epitaxial  layer  and  ion  implanted  n+-p  HgCdTe  diodes.  They 
obtained  good  agreement  with  the  experimental  data  in  general,  but  were  a 
factor  20  below  the  measured  values  at  positive  gate  bias,  when  an 
inversion  layer  formed  at  the  surface.  This  discrepancy  may  be  due  both 
to  the  presence  at  positive  gate  bias  of  a  noisy  surface  GR  contribution 
(Eqs.  120-121),  and  to  kinetic  energies  of  the  tunneling  carriers  above  the 
thermal  level  in  the  vicinity  of  the  inversion  layer. 

Another  previously  unexplained  fact  noted  was  the  difference  in  the 
fractional  noise  level  of  surface  and  bulk  recombination  currents.  This  is 
caused  in  Eqs.  (120-121)  both  by  the  surface  potential  jump  2U  of  the 
order  of  1  Volt  present  at  the  interface  between  the  bulk  and  the  oxide  and 
passivation  layers,  and  by  the  quantum  1/f  mobility  fluctuation  noise  in 
the  spreading  resistance  which  affects  the  passage  of  carriers  to  and 
from  the  perimeter  of  the  junction.  Furthermore,  the  higher  noise  level  of 
ZnS  -  passivated  diodes  may  be  caused  by  a  larger  surface  recombination 
speed  associated  with  these  coatings  compared  to  SiC>2  passivations,  and 

by  a  larger  effective  value  of  U.  The  larger  surface  recombination  speed 
pulls  more  of  the  recombination  current  from  the  bulk  to  the  surface 
where  it  has  higher  fractional  noise.  The  larger  potential  jump  U 
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increases  the  applicable  Hooge  parameter  according  to  Eq.  (121).  Finally, 
the  larger  fractional  1/f  noise  levels  of  ion  implanted  junctions  is  mainly 
caused  by  the  1-2  orders  of  magnitude  lower  carrier  lifetimes  in  Eqs. 

(1 1 8)-(1 32),  which  yield  1-2  orders  of  magnitude  smaller  Neff  values  and 

larger  fractional  noise  power  values  by  the  same  factor. 

In  order  to  reduce  the  fractional  noise  level,  our  theory  suggests  the 
use  of  a  surface  passivation  which  lowers  the  surface  recombination 
speed  and  the  surface  potential  jump  U.  The  ideal  "surface"  would  be  a 
gradual  increase  of  the  gap  width  starting  from  the  bulk  through 
compositional  changes  leading  to  a  completely  insulating  stable  surface 
outwards,  without  the  generation  of  surface  recombination  centers.  In 
addition,  the  life  time  of  the  carriers  should  be  kept  high,  and  abrupt  or 
pinched  regions  in  the  junction  should  be  avoided.  The  reasonable  choice 
of  other  junction  parameters,  including  the  steepness  of  the  junction  and 
the  geometry  should  yield  lower  injection  -  extraction  and  bulk 
recombination  noises  by  emphasizing  the  presence  of  the  larger  hole 
masses  in  the  denominators  of  the  above  expressions.  Finally,  coherent 
state  1/f  noise  should  be  avoided  by  all  means  by  optimizing  the 
dimensions,  i.e.,  keeping  them  sufficiently  small  for  individual  devices. 
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111.4  QUANTUM  1/f  NOISE  IN  MIS  INFRARED  DETECTORS 

A.  INTRODUCTION 

The  quantum  1/f  noise  theory  is  applied  here  to  the  calculation  of  1/f 
noise  in  the  various  dark  current  components  which  limit  the  performance 
of  MIS  infrared  detectors;  the  photo-current  carries  no  quantum  1/f  noise, 
because  the  elementary  process  of  photo-generation  of  an  electron  -  hole 
pair  is  almost  free  of  infrared  radiative  corrections,  and  later-on  the 
carriers  are  collected  anyway.  The  dark  diffusion  current  transports 
minority  carriers  from  the  bulk  to  the  surface  inversion  layer  and  is 
affected  by  coherent  states  quantum  1/f  noise  in  the  mobility,  as  long  as 
the  device  is  larger  than  approximately  10  pm2  in  area.  For  smaller 
devices  we  expect  the  smaller  conventional  quantum  1/f  noise  of  the 
scattering  cross  sections  to  be  expressed  in  the  mobility  and  diffusion 
fluctuations.  Quantum  1/f  recombination  speed  fluctuations  appear  as 
fluctuations  in  the  thermally  generated  current  of  minority  carriers  both 
at  the  surface  and  in  the  depletion  region,  which  are  majority  carriers  in 
the  inversion  layer.  The  rate  of  tunneling  also  presents  quantum  1/f 
fluctuations  which  are  calculated  both  for  band  to  band  tunneling  and  for 
tunneling  via  ba.'.dgap  states,  with  the  help  of  the  same  fundamental 
quantum  1/f  formula  used  for  the  diffusion  and  recombination  currents. 
Conventional  quantum  1/f  contributions  are  smaller  for  holes  than  for 
electrons,  because  they  are  inversely  proportional  to  the  effective  mass. 

B.  CURRENTS  IN  MIS  DETECTOR  STRUCTURES 

MIS  detectors^  are  different  from  photovoltaic  (or  pn  junction) 
detectors,  because  they  do  not  contain  a  pn  junction  obtained  by 
inhomogeneous  doping,  and  use  an  insulated  field  plate,  or  gate,  placed  on 
top  of  a  homogeneously  doped  narrow-bandgap  semiconductor.  The  gate  is 
used  to  control  the  surface  potential,  driving  the  semiconductor  surface 
into  deep  inversion.  The  field  of  the  induced  quasi-pn  junction  obtained 
under  the  surface  of  the  gate  in  the  homogeneous  semiconductor  material 
is  used  to  separate  the  carriers  generated  by  photo-electrically  induced 
band  to  band  transitions  just  as  in  a  photovoltaic  device.  The  MIS  infrared 
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detector  is  therefore  similar  to  a  capacitively  coupled  photovoltaic 
detector,  without  the  inconvenience  of  inhomogeneous  doping  processes. 

MIS  detectors  are  operated  in  the  pulsed  regime  by  applying  the  gate 
potential  which  creates  the  inversion  under  the  surface  for  a  finite  time 
only,  and  by  applying  subsequently  a  potential  which  flattens  the  energy 
bands  near  the  surface  and  releases  the  carriers  which  had  accumulated 
from  photo-electric  effect  and  dark  current  processes  during  the 
preceding  interval.  The  electrical  signal  obtained  when  the  carriers  are 
released,  i.e.,  during  readout,  is  proportional  to  the  number  of  carriers 
accumulated,  and  therefore  to  the  total  current  supplying  the  inverted 
volume  under  the  surface  with  minority  carriers  from  the  bulk  and  from 
various  thermal  and  photo-electric  processes  in  the  depletion,  inversion 
and  surface  regions.  This  electrical  signal  is  used  in  order  to  determine 
the  flux  of  infrared  radiation.  For  this  determination,  however,  the  dark 
current  contribution  needs  to  be  subtracted  first. 

The  dark  current  is  the  current  supplying  the  potential  well,  defined  by 
the  inversion  region  under  the  surface,  with  minority  carriers  in  the 
absence  of  the  applied  infrared  flux.  Any  low  -  frequency  fluctuation  in 
the  dark  current  will  be  interpreted  as  a  fluctuation  in  the  major  infrared 
flux  signal.  Therefore  fluctuations  of  the  dark  current  at  frequencies 
below  the  readout  frequency  will  limit  the  performance  of  infrared 
detectors.  In  the  pulsed  mode  of  operation  considered  here,  the  dark 
current  is  monitored  only  during  the  inverted  phase,  when  carriers  are 
accumulated  in  the  potential  well.  Therefore  the  cross  sections  and 
process  rates  which  control  the  intensity  of  the  dark  current  are  not 
sampled  continuously  by  us  either.  Nevertheless,  the  quantum  1/f 
fluctuations  of  these  cross  sections  and  process  rates  will  be  practically 
the  same  as  if  we  would  have  observed  them  continuously.  Indeed,  we 
first  note  that  the  changes  caused  by  the  applied  gate  voltages  in  the 
microscopic  flux  of  electrons  testing  all  cross  sections  and  process  rates 
defined  in  a  semiconductor  do  not  affect  the  cross  sections  or  process 
rates  themselves.  There  is  a  similar  flux  interrogating  the  cross  sections 
in  thermal  equilibrium.  How  much  and  in  which  way  this  flux  is  modified 
by  the  applied  voltage  is  not  important;  there  is  always  some  flux  there  to 
test  the  cross  sections  at  all  energies  and  scattering  angles.  Being  the 
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ratio  of  the  scattered  current  to  the  incoming  flux,  the  cross  section  is 
not  affected  by  flux  changes.  Second,  we  note  that  since  it  contains  the 
number  of  carriers  N  in  the  denominator,  quantum  1/f  noise  should  be 
dependent  in  general  on  the  flux  level  used  in  testing  the  cross  sections  or 
process  rates.  However,  N  is  not  the  number  of  carriers  accumulated  in 
the  well  during  the  pulse  of  the  gate  voltage.  N  is  the  number  of  carriers 
used  to  define  the  scattered,  or  outgoing  current  at  any  time  during  the 
pulse  of  the  gate  voltage,  i.e.,  the  number  of  carriers  effectively  carrying 
the  scattered  current  in  stationary  scattering  conditions.  The  stationary 
scattering  conditions  are  established  in  a  time  of  the  order  of  the  mean 
carrier  collision  time,  which  is  much  shorter  than  the  pulse  length  of  the 
gate  voltage.  Therefore,  the  value  of  N  during  pulsed  operation  is  not 
dependent  on  the  pulse  length,  and  is  not  different  from  the  value 
corresponding  to  an  infinite  pulse  length  or  continuous  (cw)  operation, 
even  if  the  applied  voltage  is  so  large  that  the  microscopic  flux 
interrogating  the  cross  sections  differs  considerably  from  the  thermal 
equilibrium  value.  We  conclude  that,  indeed,  the  quantum  1/f  noise 
affecting  the  cross  sections  and  process  rates  is  not  changed  by  the 
pulsed  operation.  This  independence  of  1/f  noise  on  the  continuous  or 
discontinuous  character  of  any  applied  bias  has  been  experimentally 
verified38  during  the  last  2  decades,  and  has  been  found  to  be  in 
agreement  with  the  interpretation  of  1/f  noise  in  terms  of  fundamental 
resistance  fluctuations.  Although  the  experimental  verification  was 
performed  on  fluctuations  in  conduction  only  (involving  mainly  scattering 
cross  sections),  from  the  above  discussion  based  on  the  concept  of 
quantum  1/f  noise  we  know  that  the  similarity  of  quantum  1/f  noise  in 
the  continuous  and  pulsed  regimes  should  be  also  true  for  quantum  1/f 
fluctuations  in  recombination  cross  sections  and  tunneling  rates. 

The  dark  current  has  to  be  subtracted  from  the  total  current  in  a 
(HgCd)Te  MIS  device  to  yield  the  photocurrent.  Therefore,  the  minority  - 
carrier  dark  current  is  the  single  most  important  parameter  for  the 
operation  of  MIS  devices  as  detectors  of  infrared  radiation33.  This 
applies  both  to  operation  of  MIS  devices  in  the  thermal  equilibrium  mode, 
in  which  the  dark  current  determines  the  MIS  diode  impedance,  and  to 
operation  in  the  dynamic,  or  integrating  mode,  in  which  the  gate  voltage  is 
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pulsed,  and  in  which  the  minority  -  carrier  dark  current  determines  the 
storage  time  of  the  device.  The  main  component  of  the  dark  current  in 
narrow-bandgap  HgCdTe  is  the  tunneling  current  via  bandgap  states35, 
which  can  also  be  considered  as  an  electric  breakdown  effect.  In  general, 
the  tunneling  current  occurs  both  through  band  to  band  transitions  and 
through  intermediary  states.  The  band  to  band  tunneling  current  through  a 
simple  trianguiar  barrier  is 

ltb  -  (q3EOs/4ji3h2)(2mVEg)1/2exp[4(2m*)^2Eg3/2/3qbE],  (134) 

where  E  is  the  electric  field  associated  with  the  barrier,  and  Eg  is  the 
band  gap.  The  electric  field  can  be  approximated  by  the  electric  field  at 
the  semiconductor  surface 

E  -  (2qn0<I>s/e£o)1/2.  (135) 

where  <DS  represents  the  empty  well  surface  potential,  and  n0  is  the  doping 
concentration.  Substituting  this  value  into  Eq.  (135),  with  m*/m0  *7  10' 2 

Eg,  we  obtain3  5 

ltb  -  10-2no1/2<Ds3/2exp[-  4.3  10l0Eg2/(no<I>s)1/2]  A/cm2,  (136) 

where  n0  is  in  crrr3  and  Eg,  <DS  in  volts.  Therefore  the  tunneling  current  is 
strongly  dependent  on  the  bandgap  and  also  depends  on  the  doping 
concentration  and  the  surface  potential. 

Experimental  values  of  the  tunneling  current  are  usually  larger  than  Eq. 
(136)  because  of  the  additional  effect  of  tunneling  via  bandgap  states. 

This  effect  is  particularly  important  in  n  -  type  devices39.  Indeed,  in  n  - 
type  devices  the  applied  gate  voltage  is  negative  in  order  to  produce 
depletion  at  the  surface.  The  energy  bands  are  therefore  curved  upwards 
at  the  surface,  and  transitions  of  electrons  from  the  valence  band  to 
Shockley  -  Read  (SR)  states  at  the  middle  of  the  band  gap,  as  well  as  the 
subsequent  transitions  from  these  states  to  the  conduction  band  are 
facilitated  by  the  presence  of  many  defects  right  at  the  surface  of  the 
semiconductor.  In  p  -  type  devices  the  similar  indirect  tunneling 
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processes  occur  farther  away  from  the  surface,  because  in  this  case  the 
bands  are  curved  downwards  at  the  surface,  and  transitions  of  electrons 
from  the  valence  band  to  the  centers  at  the  middle  of  the  band  gap,  as  well 
as  the  transitions  from  the  centers  to  the  conduction  band  well  at  the 
surface,  occur  right  where  the  curvature  begins,  i.e.,  further  away  from 
the  surface.  We  conclude  that  in  p  -  type  devices  there  will  be  fewer  SR 
centers  active  in  indirect  tunneling,  and  therefore  the  tunneling  current  lc 

via  SR  centers  at  a  given  temperature  and  a  given  applied  gate  voltage 
will  be  smaller.  The  tunneling  current  will  be  further  reduced  in  p  -  type 
devices  due  to  the  lower  density  of  states  present  in  the  surface  potential 
well  due  to  quantization  of  the  motion  of  the  electrons  in  the  potential 
well  at  the  surface.  The  reduced  values  of  the  dark  current  in  p  -  type 
devices  correspond  to  higher  values  of  the  breakdown  field  in  these 
devices.  The  best  measured  value39  of  the  breakdown  field  in  10  pm 
cutoff  p  -  type  devices  is  in  excess  of  1.0  V/m,  whereas  that  for  n  -  type 
material  of  similar  bulk  defect  quality  is  0.5  V/m.  On  the  other  hand,  the 
minority  carriers  diffusion  current  is  larger  in  p  -  type  devices  due  to  the 
smaller  mass  and  higher  diffusion  constant  of  electrons  compared  to 
holes.  The  advantage  of  p  -  type  devices  is  therefore  considerable  only  in 
the  case  of  very  narrow  band  gap  and  very  long  cutoff  wavelengths.  We 
shall  therefore  consider  both  the  case  of  p  -  type  and  n  -  type  devices. 

The  large  diffusion  current  present  in  p  -  type  devices  corresponds  to  the 
large  value  of  the  diffusion  length  of  electrons  and  can  be  reduced  by 
thinning  the  device,  i.e.,  by  reducing  its  thickness  well  below  the  diffusion 
length. 

In  general  the  dark  current  I  can  be  written  in  the  form 

I  »  Id  +  lr  +  Is  +  Ifb  +  Itc  +  Ifsc  +  lb  +  cpl®-  0  37) 

The  eight  terms  on  the  right  hand  side  correspond  to  minority  carrier 
diffusion  from  the  bulk,  generation-recombination  from  SR  centers  in  the 
depletion  region,  generation  from  SR  centers  at  the  surface,  band  to  band 
tunneling,  tunneling  via  SR  centers,  tunneling  via  surface  centers, 
recombination  on  the  back  surface  and  photoelectric  generation  with 
efficiency  r\  by  the  thermal  radiation  background  flux  O.  We  shall  give  the 
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formulae39  which  determine  these  terms  below,  also  including  an  example 
of  their  calculation  in  a  p  -  type  device 

ld  -  (qnj2/p0)[kWqxn]1/2  -  0-6  10‘19c  36  1024  cm-6/1 0*1 5cm-3) 
[0.01V  1.5  105(cm2/Vs)/10-6s]1/2  .  2  10'4A/cm2,  (138) 

lr  -  qnjW/2x0  -  [(1.6  lO'^C  6  1012cm-3  2  10-4cm  /(2  lO'^s)] 

-  10-4A/cm2  (139) 

ls  (=  lb  )=  qnjS/2  -  [(1.6  10‘19c  6  1012cm-3  20cm/s] 

-  1.8  10‘5A/cm2,  (140) 

ltc  .  lO-^NrOs/EgJexpt-S.S  106Eg2/E]  -10-'|3(4  1012cm-3o.2V/0.062V) 
exp[-5.3  106(0.062V)2/3000V/cm]  -  8  10'4  A/cm2,  (141) 

and  Itb  was  given  by  Eq.  (136),  yielding  10‘7A/cm2.  Here  nj  -  6  1012cm3 

is  the  intrinsic  carrier  concentration,  and  the  concentration  of  holes  was 
taken  to  be  p0  -  1015cm_3.  The  mobility  p  *  1.5  1015cm2/Vs  as  well  as 
the  life  time  xn  »  10_6s  of  the  minority  carriers  have  to  be  replaced  by  pp 
and  tp  for  the  case  of  n  -  type  devices.  The  surface  recombination  speed 
was  denoted  by  s  -  40cm/s,  and  the  concentration  of  intermediary  states 
effective  for  tunneling  was  denoted  by  Nr  =  4  10l2cm-3.  The  band  gap 
considered  was  Eg  -  0.065eV,  the  surface  potential  Os  »  0.2V  and  the 
electric  field  below  the  surface  E  -  3000V/cm.  All  numerical  values  have 
been  included  only  as  an  example  and  are  not  characteristic  of  a  particular 
device.  The  numerical  factors  in  Eq.  (141)  correspond  to  p  -  type  HgCdTe 
with  a  20  pm  cutoff  wavelength  and  were  taken  from  the  paper  of  Kinch 

and  Beck39.  For  n  -  type  devices  we  also  need  to  include  tunneling  via 
surface  states.  If  the  density  of  fast  surface  states  is  denoted  by  Nfs  = 

1012  cm*2v1,  the  current  generated  from  tunneling  via  a  uniform 
density  of  fast  surface  states  across  the  band  gap  will  be  given  by35 


Itsc  *  -qPlNfskT/2q  •  1.75  10-4  A/cm2. 


(142) 
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The  sum  of  the  first  seven  currents  on  the  right  hand  side  of  Eq.  (137) 
must  be  smaller  than  the  eighth  which  corresponds  to  the  thermal 
radiation  background  flux,  for  background  limited  (BLIP)  performance. 
Although  not  all  terms  in  the  dark  current  are  of  importance,  we  still 
retain  them  at  this  point,  because  their  quantum  1/f  noise  may  be  quite 
significant,  even  if  the  corresponding  current  is  negligible.  We  shall  now 
proceed  with  the  calculation  of  quantum  1/f  noise  contributions  from  all 
these  currents. 

C.  QUANTUM  1/f  NOISE  SOURCES 

1.  1/f  Noise  in  the  Diffusion  Current 

The  diffusion  limited  dark  current  l<j  will  exhibit  1/f  noise  due  to 
conventional  quantum  1/f  fluctuations  in  the  scattering  cross  sections  of 
the  carriers  due  to  phonons  and  impurities.  We  apply  the  fundamental 
formula  given  by  Eq.  (3)  for  an  individual  scattering  process  in  which  the 
velocity  change  v  is  given  by  the  thermal  energy  of  the  carriers,  with  the 
assumption  that  the  collisions  are  perfectly  randomizing  collisions.  If 
the  velocity  v  is  rotated  by  an  angle  8  in  an  elastic  collision,  the  velocity 
change  magnitude  is  |Av|  »  2v  sin(0/2).  Averaging  over  all  angles  and 
velocities,  we  obtain 

Av2  =  4v2  sin2(e/2)  *  2v2,  (143) 

and  therefore  from  Eq.  (1.1)  we  get  in  thermal  equilibrium  at  the 
temperature  T  the  1/f  noise  coefficient 

aH  -  (4a/37t)(6kT/m‘c2),  (144) 

where  we  have  assumed  a  Maxwell  distribution  of  velocities.  For  Hg-j.x 
CdxTe  with  x  -  0.2  we  have  mn*  *  0.008mo  and  for  x  *  0.3  we  have  mn*  = 
0.02mo.  Therefore  we  obtain  an  -  2  10"7  in  the  first  case  and  an  -  7.5 
10-8  in  the  second  case. 
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For  the  case  of  Umklapp  scattering,  which  occurs  in  semiconductors 
only  to  a  limited  extent  due  to  the  relatively  small  number  of  high 
momentum  phonons  available  at  the  temperature  T,  the  momentum  change 
of  the  electron  is  given  by  the  smallest  reciprocal  lattice  vector,  and 
therefore  |Av|  =  h/am*.  We  therefore  obtain  the  quantum  1/f  noise 
coefficient 

an  *  (4a/37c)(h/m*ac)2,  (145) 

which  is  much  larger  than  Eq.  (144),  but  has  to  be  multiplied  with  a 
negative  exponential  which  describes  the  scarcity  of  phonons  with 
momentum  of  the  order  of  a  reciprocal  lattice  vector.  [The  negative 
exponential  e*®/T  could  be  included  in  the  current  weight  factor  which 
will  be  defined  below  in  Eq.  (),  but  we  prefer  to  include  it  here  already]. 
Combining  Eqs.  (144)  and  (145),  we  obtain  for  conventional  1/f  noise  in 
the  mobility  and  diffusion  coefficients 

aH  =  (4a/3ic)[(6kT/m*c2)  +  (h/m*ac)2exp(-0/T)],  (146) 

where  0  is  about  half  the  Debye  temperature  for  simple  metals,  but  may 
be  higher,  of  the  order  of  the  Debye  temperature,  for  semiconductors. 

The  quantum  1/f  noise  considered  so  far  is  known  as  conventional 
quantum  1/f  noise,  and  affects  cross  sections  and  process  rates.  In 
sufficiently  large  semiconductors  samples  we  expect  coherent  state 
quantum  1/f  noise.  For  this  type  the  1/f  noise  coefficient  is  given  by 

«coh  -  2a/jt  =  4.6  10'3.  (147) 

The  values  of  the  quantum  1/f  noise  coefficient  given  by  Eqs.  (143)  - 
(147)  can  be  used  to  calculate  the  quantum  1/f  noise  which  affects  the 
various  currents  listed  in  Eq.  (137).  We  first  consider  the  case  of  the  dark 
diffusion  current  of  electrons  from  the  bulk  through  the  surface  barrier  in 
a  p  -type  MIS  device,  similar  to  diffusion  in  a  n+p  junction,  because  in 
both  cases  the  current  is  determined  by  the  diffusion  of  electrons  which 
are  minority  carriers,  against  the  built  -  in  field  of  a  Boltzmann  potential 
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barrier  into  the  surface  well,  and  by  the  thermal  generation  of  carriers 
there.  We  start  with  the  derivation  of  the  mobility  fluctuation  part  of 
quantum  1/f  noise  in  a  n+-p  diode.  For  the  MIS  barrier,  just  as  for  a 
diffusion  limited  n+-p  junction,  the  current  is  controlled  by  diffusion  of 
electrons  into  the  p  -  region  over  a  distance  of  the  order  of  the  diffusion 
length  L  -  (Dnxn)1/2  which  is  usually  shorter  than  the  length  wp  of  the  p  - 
region.  If  N(x)  is  the  number  of  electrons  per  unit  length  and  Dn  their 
diffusion  constant,  the  electron  current  at  x  is 


lnd  *  *  eDndN/dx, 


(148) 


where  we  have  assumed  a  planar  junction  and  taken  the  origin  x  =  0  in  the 
junction  plane.  Diffusion  constant  fluctuations,  given  by  kT/e  times  the 
mobility  fluctuations,  will  lead  to  local  current  fluctuations  in  the 
interval  Ax 

5Alnd(x,t)  «  lnd^x8Dn(x,t)/Dn.  (149) 

The  normalized  weight  with  which  these  local  fluctuations  representative 
of  the  interval  Ax  contribute  to  the  total  current  Id  through  the  diode  at  x 

«  0  is  determined  by  the  appropriate  Green  function  and  can  be  shown  to  be 
(1/L)exp(-x/L)  for  wp/L  »  1.  Therefore  the  contribution  of  the  section  Ax 

i  s 


5Ald(x,t)  =  (Ax/L)exp(-x/L)lnd8Dn(x,t)/Dn,  (150) 

with  the  spectral  density 

Sa|d(x.f)  -  (Ax/L)2  exp(-2x/L)  lnd2  SDr,(x,f)/Dn2.  (151) 

For  mobility  and  diffusion  fluctuations  the  fractional  spectral  density  is 
given  by  annd^Nx.  where  annd  IS  determined  from  quantum  1/f  theory 
according  to  Eqs.  (143)  -(147).  With  Eq.  (148)  we  obtain  then 


SA|d(x,f)  -  (Ax/L2)  exp(-2x/L)  (eDndN/dx)2  aHnd/fN- 


(152) 
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The  electrons  are  distributed  according  to  the  solution  of  the  diffusion 
equation,  i.e. 

N(x)  -  [N(0)  -Np]  exp(-x/L)  +Np;  dN/dx  -  -{[N(0)  -Np]/L>  exp(-x/L).  (153) 

Substituting  into  Eq.  (152)  and  simply  summing  over  the  uncorrelated 
contributions  of  all  intervals  Ax,  we  obtain 

wp 

Sld(f)  -  aHnd(eDn/L2)2  J[N{0)  -  Np]2e-4x/Ldx  /{[N(0)  -  Np]e-*'L  +Np}.  (154) 

0 

We  note  that  (eDn/L2)  *  (e/xn).  With  the  expression  of  the  saturation 
current  lo  -  e(Dn/n)1/2Np  and  of  the  current  I  -  lo[exp(eV/kT)  -  1]  we  can 
carry  out  the  integration 

1 

S|d(r)  -  aHnd^^n)  Ja2u3du/(au  +  1) 

0 

-  aHnd(el/ftn)  [F(a)  -F(a  e-wp/L)]  *  aHnd(el/ftn)awp/[(a  +  1)L],  (155) 

the  last  form  being  for  wp  «  L.  Here  F(a),  a  and  u  are  given  by  Eq.  (93). 

For  wp  »  L  we  have  F(o)  -  0,  and  the  second  term  in  rectangular  brackets 
drops  out  in  Eq.  (155). 

Eq.  (155)  only  contains  the  fluctuations  in  the  mobility  and  the 
diffusion  constant.  In  a  similar  way  we  calculate  the  quantum  1/f 
fluctuations  of  the  recombination  rate  in  the  bulk  of  the  p  -  region.  We 
have  for  the  recombination  current  AIr(x)  in  a  section  Ax,  if  N'(x)  is  the 

excess  carrier  density, 

AIr(x)  -  eN'(x)Ax/xn  (156) 

Putting  Cn  -  1/xn  and  bearing  in  mind  that  xn,  and  hence  Cn,  fluctuates,  we 
have  for  the  section  Ax, 

6AIr(x,  t)  =  AlR(x)[5Cn/Cn]  (157) 

with 
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SC„(x,f)/Cn2  -  ahnr/fN,  (158) 

so  that,  since  N(x)  ®  N/Ax, 

S4|r(x,  f)  »  AIR2(x)[aHn,/fN]  -  aHnr(e2[N'(x)]2/ftn2N(x)}Ax  (159) 

where  N'(x)  »  [N(0)  -  Np]exp(-x/L)  and  N(x)  »  N’(x)  +  Np  as  before. 

It  is  easily  shown  that  the  fluctuating  current  Al(x,  t)  at  the  junction 
i  s 


Al(x,  t)  =  AIr(x,  t)exp(-x/L),  (160) 

so 

wp 

S,(f)  -  aHnr(e2NpL/ftn2)  J{[N'<x)j2/NpN(x)}exp(-2x/L)d(x/L) 

0 

1 

=  aHnr(el0/fxn)  J(a2u3/(au  +  1)]du  -  oH„r[el/ftn][F(a)  -  F(ae-»P'L)],  (161) 
0 

where  F(a),  a,  and  u  have  the  same  meaning  as  before. 

We  can  use  the  similarity  of  the  quantum  1/f  noise  results  for 
diffusion  current  fluctuations  caused  by  mobility  fluctuations  and  by 
recombination  speed  fluctuations  in  order  to  combine  both  into  a  single 
formula 

S|(0  -  (aHd  +  «Hnr)[el/ftn](F(a)  -  F(ae-wp/L)].  (162) 

In  the  limit  of  very  short  devices  (wp  «  L)  the  last  factor  becomes 
awp/[(a  +  1)L],  and  in  the  limit  of  long  MIS  devices  (wp  »  L)  it  simply 
becomes  F(a).  In  addition  we  have  a  current  noise  contribution  S|b  from 
the  quantum  1/f  fluctuation  of  the  recombination  speed  s  on  the  back 
surface. 

So  far  we  have  const  Ted  only  conventional  quantum  1/f  noise  which 
is  applicable  to  sufficiently  small  devices.  In  general,  however,  we  must 
interpolate  between  conventional  and  coherent  quantum  1/f  noise, 
according  to  the  relation 
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ccH  -  [1/(1  +  s)][2aA/f]  +  [s/(1  +  s)][2a/rcf],  (163) 

where  s  -  Em  /E|<  -  2e2N7m*c2,  and  N’  -  nS  is  the  number  of  the  carriers 
per  unit  length  of  the  device.  Here  s  represents  the  ratio  between  the 
magnetic  energy  per  unit  length  and  the  kinetic  energy  per  unit  length  of 
the  device.  The  quantity  e2/m*c2  -  r0m/m*  can  be  calculated  in  terms  of 
the  classical  radius  of  the  electron  r0  «  2.84  10'13  cm.  Then  we  obtain  s 
-  2  r0N'm/m\  i.e.,  the  parameter  s  represents  twice  the  number  of 
carriers  present  in  a  length  of  the  device  equal  to  the  classical  radius  of 
the  electron.  We  must  compare  s  with  A,  and  if  s  «  A  we  apply 
conventional  quantum  1/f  noise,  whereas  for  s  »  A  we  have  to  apply 
coherent  quantum  1/f  noise.  In  genera!  the  approximate  formula  of  Eq. 
(163)  must  be  used  for  the  transition  region. 

The  dimensionless  parameter  s  is  easy  to  calculate  in  any  practical 
case.  For  instance  in  the  case  of  a  MIS  device  of  area  50  pm  x  50  pm  with 
a  concentration  of  carriers  of  1016  cm-3  we  obtain  N'  =  2.5  lO^/cm,  and 
with  m/m*  =  50  we  obtain  s  =  7.  On  the  other  hand,  we  can  estimate  A  for 
conventional  quantum  1/f  noise  and  we  will  certainly  find  A  «  1,  because 
the  velocity  change  of  the  carriers  must  be  much  smaller  than  the  speed 
of  light.  Therefore,  in  this  case  we  must  apply  coherent  quantum  1/f 
noise,  because  s  »  A.  Consequently,  in  Eq.  (3.20)  we  must  set 

«Hd  +  «Hnr  -  «coh  -  4.6  10'3  (164) 

The  coherent  state  quantum  1/f  noise  coefficient  thus  replaces  the  total 
conventional  Hooge  parameter. 

2.  1/f  Noise  of  the  Recombination  Current  Generated  in  the 

Depletion  Region 

The  quantum  1/f  noise  of  the  recombination  current  thermally 
generate  jn  the  depletion  region  arises  from  quantum  1/f  fluctuations  of 
the  bulk  recombination  rates  in  the  depletion  region.  The  difference 
between  the  recombination  rate  R  and  the  generation  rate  G  is  given  by 
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R(x)  -  G(x)  -  [pn  -  n2j]/[(n  +  n-|)xp0  +  (p  +  p-|)xno]  (165) 

where  ni  and  p-j  are  electron  and  hole  densities  when  the  Fermi  level  lies 
at  the  trap  level.  If  the  trap  level  lies  at  the  intrinsic  level,  n-j  *  p-j  *  nj. 
Moreover,  xpo  and  xno  are  time  constants  for  electrons  and  holes.  If  A  is 
the  cross  -  sectional  area  of  the  junction,  the  current  is 

w  w 

lr  -  e  J[R(x)  -  G(x)]A  dx  =  e  J[pn  -  n2j3/[(n  +  nj)xpo  +  (p  +  nj)xno]A  dx,  (166) 
0  0 

where  w  is  the  width  of  the  space  -  charge  region  and  the  trap  level  is 
assumed  to  lie  at  the  intrinsic  level. 

We  now  turn  to  the  g  -  r  noise.  The  time  constants  xpo  and  xno  fluctuate 
in  a  1/f  fashion  due  to  the  quantum  1/f  fluctuations  in  the  recombination 
cross  sections,  and  this  produces  the  quantum  1/f  contribution  to  g  -  r 
noise.  We  now  write 

xno  *  1/C  ni  xpo  *  1/Cp,  (167) 

where  Cn  and  Cp  are  the  generation  (or  combination)  rates  for  a  single 
electron  and  for  a  single  hole,  respectively.  Consequently 

SxncAno  ■  -(5Cn/Cn);  5xpo/xpo  *  *(5Cp/Cp)  068) 

We  now  apply  this  knowledge  to  Eq.  (165)  and  observe  that 

8(R  -  G)  *  [R(x)  -  G(x)]x 

{[(n  +nj)xpo(8Cp/Cp)]  +[(p  +nj)xno(8Cn/Cn)]}/[(n  +nj)xpo  +  (p  +nj)xno].  (169) 
so  that  with 


w 

8lr  «  e  Js[R(x)  -  G(x)]Adx, 
0 


(170) 
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the  noise  is 

w  w 

S|,(0  -  e2  }  /  (R(x)  -  G(x)]A[(R(x')  -  G(x’)]A 
0  0 

{[(n  +  nj)2tpo2Scp(x,  x',  f)/Cp2  +  (p  +  nj)2xno2Scn(x.  x\  f)/Cn2] 

/[(n  +  nj)tpo  +  (p  +  nj)tno]2}dxdx\  (171) 

since  Cp  and  Cn  are  independent. 

We  now  observe  that 

pn  -  n2j  -  n2j[exp(eV/2kT)  -  1][exp(eV/2kT)  +  1],  (172) 

and  that  the  integrand  in  Eq.  (171)  has  an  appreciable  value  only  if  p  «  n  * 
nj  exp(eV/2kT).  By  substituting  n  +  nj  =  p  +  nj  *  nj[exp(eV/2kT)  +  1]  we 
define  an  effective  width  weff  such  that 

w 

eA  J[(pn  -  n2j)dx]/[(n  +  nj)xpo  +  (p  +  nj)tno] 

0 

-  eA{(pn  -  n2j)/nj[exp(eV/2kT)  +  1]}[weff/(xpo  +  tno)l-  (173) 
We  may  thus  write 

lr  -  lgr  -  eAweff n j[exp(e V/2kT)  -  1]/(xpo  +  xno) 

-  [eNeff/(Tpo  +  xno)Jtanh  eV/2kT,  (174) 

where  Neff  »  Aweffnj[exp(eV/2kT)  +  1]  is  the  effective  number  of  hole  - 
electron  pairs  taking  part  in  the  conduction  and  noise  processes.  This 
equation  is  exact  but  not  very  useful  since  it  contains  the  unknown 
parameter  weff. 

We  now  turn  to  Eq.  (171)  and  observe  that 


Scp(x.  X',  f)/C2p  -  {(aHp/f)/[R(x')  +  G(x')](Tp0  +  tn0)A}8(x'  -  x),  (175) 
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and 

SCn(x,  x*,  f)/C2n  -  {(aHn/f)/[R(x')  +  G(x')](xpo  +  tno)A}5(x*  -  x).  (1 76) 

The  factor  (xno  +  xpo)  comes  in  because  Scp(x,  x\  f)/Cp2  and  Scn(x,  x‘, 
f)/Cn2  must  be  independent  of  xp0  and  xno  if  p*n*nj  exp(eV/2kT).  This 
yields,  if  we  integrate  over  the  8  function 

w 

S|r(f)  -  e2  J{[R(x)  -  G(x)]2A/(tpo  +  ^nc)tR(x)  +  G(x)]} 

0 

{[(n  +  nj)2xpo2aHp/f  +  (p  +  nj)2xno2aHn/f]/[(n  +  nj)xpo  +  (p  +  nj)xno]2}dx  (177) 

We  now  observe  that  the  second  factor  in  Eq.  (177)  is  practically  a 
constant  as  long  as  p  and  n  are  comparable.  We  may  thus  bring  that  factor 
outside  the  integral  sign  and  write 

[(n  +nj)2xpo2aHp/f  +  (p  +nj)2xno2aHn/f]/[(n  +nj)xpo  +  (p  +n-,)xno]2  *  aH/f,  (178) 
where  an  is  given  by 

aH  “  frpo^po  +  Tno)]2<xHp  +  [^noA^po  +  tno)]2(XHn  (1  79) 


We  thus  have 

w 

S|r(f)  -  [aHe/f(tn0  +  Tp0)]e  J[R(x)  -  G(x)]2A/[R(x)  +  G(x)]dx 

0 

-  [aHelgr/f(Tno  +  XpoWanh  [eV/2kT).  ( 1 80) 

We  can  now  prove  Eq.  (180)  in  a  Hooge  -  type  formulation.  Here  we  put 
S|r(f)/Ir2  -  aH/fNeff  (181) 

But  according  to  Eq.  (174),  lgr  -  (eNeff/x)tanh  eV/2kT,  so  that 

S|(f)  -  aH[elr/fx]tanh  eV/2kT,  x  -  (xno  +  xpo)  (182) 
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in  agreement  with  Eq.  (180). 

We  can  also  prove  Eq.  (180)  from  the  following  consideration.  We 
write  lr  -  (Neffl')  tanh  eV/2kT,  where  I’  ■  e/t  fluctuates.  We  then  have 

S|'(f)/I’2  -  aH/f  or  S|-(f)  -  aH/f(e/x)l  (183) 

so  that,  since  the  Neff  hole  -  electron  pairs  are  independent 

S|(f)  -  NeffS p(f)  »  (eNefflVft)aH  *  [(elgr/fx)aH]tanh  eV/2kT,  (184) 

where  is  given  by  Eq.(179)  and  x  »  (xno  +  xpo). 

The  last  two  approaches  are  easily  extended  to  other  cases;  the  method 
works  as  long  as  a  time  constant  x  and  an  Neff  can  be  defined. 

We  finally  evaluate  <xhp  and  ann  from  quantum  1/f  noise  considerations: 
aHn  *  (4a/37t)(Av2n/c2)  =  (4a/3jt){[2ea(Vdif  -  V)  +  3kT]/m*nc2}>  (185) 

aHp  =  (4a/3rc)(Av2p/c2)  -  (4a/3w)|[2e(1  -  a)(Vdjf  -  V)  +  3k7]/m‘pc2|,  (186) 
and  as  a  consequence  (see  Eq.  179), 

aH  =■  (4a/3jc)[2e(Vdjf  -  V)  +  6kT]/[(m‘n)^2  +  (m*p)i/2]2c2  (187) 

The  problem  has  hereby  been  solved.  Note  that  in  Hgi-xCDxTe  with  x  =  0.3, 
m*n  -  0.02  m,  m*p  »  0.55  m,  so  that  [(m*n)1/2  +  (m*p)1/2]2  =  0.78  m,  very 
much  larger  than  m*n. 

3.  Noise  in  the  Surface  Recombination  Current 

The  surface  recombination  current  is  a  dark  current  component 
originating  from  surface  states  at  the  interface  between  the  surface 
passivation  layer  and  the  bulk.  The  recombination  cross  sections  of  these 
states  will  exhibit  quantum  1/f  noise.  The  quantum  1/f  noise  coefficient 
an  reflects  the  velocity  change  of  carriers  involved  in  this  process 
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aH  -  (4a/3jt)(2/m*c2)[3kT/2  +  eU/2  +  0.1  eV],  (1 88) 

where  U  is  the  surface  potential  jump  present  between  the  surface 
passivation  layer  and  the  bulk  even  if  no  gate  voltage  is  applied,  and  V  is 
the  applied  gate  voltage  which  we  take  with  a  coefficient  less  than  unity, 
here  for  example  with  a  weight  of  0.1. 

The  calculation  of  surface  recombination  quantum  1/f  noise  is  similar 
io  the  calculation  of  quantum  1/f  noise  from  recombination  in  the  space 
charge  region.  However,  in  this  case  the  cross  sections  are  not 
distributed  over  the  width  of  the  junction,  but  rather  are  concentrated  at 
the  surface.  Therefore  we  can  repeat  the  above  derivation  and  we  obtain 
an  expression  of  the  form 

S|s(f)  -  aHels[tanh  x]  /f(tno  +  xp0).  (x  *  eV/2kT)  (189) 

We  note  that  the  fabrication  process  of  MIS  structure  introduces  less 
bulk  defects  than  the  fabrication  of  photovoltaic  devices.  Nevertheless, 
the  fabrication  of  MIS  devices  introduces  some  defects  in  the  bulk  layer 
located  right  under  the  surface.  These  defects  will  manifest  themselves 
through  a  contribution  to  indirect  tunneling. 


4.  Quantum  1/f  Noise  in  the  Tunneling  Rate 

In  the  case  of  tunneling  from  a  surface  accumulation  layer  to  the  bulk, 
the  velocity  change  of  the  carrierc  will  lead  us  from  the  thermal  velocity 
of  the  carrier  on  one  side  of  the  barrier  to  the  thermal  velocity  of  a 
carrier  on  the  other  side  of  the  barrier,  if  band  to  band  tunneling  is 
considered.  If,  however,  the  tunneling  goes  via  intermediate  states 
located  in  the  bandgap,  the  velocity  is  zero  as  long  as  the  carrier  is 
stationary  in  the  intermediate  state.  We  can  therefore  write  the  1/f 
noise  coefficient 


an  -  (4a/37t)6kT/m*c2 


(190) 
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for  band  to  band  tunneling,  and 

<xH  -  (4a/37t)3kT/m*c2  (191) 

for  tunneling  via  intermediate  states  in  the  band  gap,  where  we  have 
considered  the  average  squared  velocity  change  two  times  smaller.  The 
effective  mass  is  the  mass  of  the  minority  carriers  in  the  bulk  material. 

For  band  to  band  tunneling  from  a  surface  inversion  layer  to  the  bulk, 
the  velocity  change  of  carriers  corresponds  to  an  energy  difference  of  the 
order  of  the  band  gap  Eg  plus  an  energy  difference  of  the  order  of  the 

thermal  energy  3kT/2,  provided  we  are  dealing  with  deep  inversion,  as 
used  in  practical  MIS  devices  in  the  pulsed  mode  of  operation.  This  yields 
the  quantum  1/f  coefficient 

<xH  -  (4a/37t)(Eg  +  3kT/2)/m‘c2  (192) 

for  band  to  band  tunneling.  For  tunneling  via  intermediate  states  in  the 
band  gap  the  corresponding  energy  difference  will  be  smaller,  and 
therefore  we  replace  Eg  by  Eg/2 

aH  -  (4a/37c)(Eg/2  +  3kT/2)m*c2  (193). 

This  relation  is  applicable  both  if  the  intermediate  states  are  located  in 
the  depletion  region  or  at  the  surface. 

In  the  last  four  equations  we  did  not  divide  by  the  number  of  carriers 
simultaneously  involved  in  the  tunneling  process,  because  this  number  is 
less  than  unity  for  practical  tunneling  currents.  Whenever  a  cross  section 
or  a  process  rate  is  tested  with  one  electron  or  less  than  one  electron  at  a 
time,  the  effective  number  of  electrons  in  the  denominator  of  the  quantum 
1/f  formula  must  be  replaced  by  unity.  Let  us  calculate  the  average 
number  of  carriers  simultaneously  present  in  the  tunneling  process  at  any 
time.  The  tunneling  process  occurs  over  a  distance  d  -  Eg/eE,  and  the 

speed  v  of  the  carriers  will  be  of  the  order  of  the  thermal  speed  in  the 
case  of  an  accumulation  layer,  and  of  the  order  of  the  band  gap  energy  in 
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the  case  of  a  surface  inversion  layer.  Dividing  d  by  v,  we  obtain  a 
tunneling  time  t  »  6  10‘12s  for  accumulation  layers  and  t  -  3  lO'^s  f0r 
inversion  layers.  From  Eq.  (141)  we  know  that  the  tunneling  current  is  of 
the  order  of  10'2A/cm2.  Multiplying  this  by  t,  we  obtain  3  -  6  lO*1®  C, 
i.e.,  2000  -  3000  electrons/cm2  tunneling  simultaneously  in  a  device  of  1 
cm2  area.  In  a  device  of  dimensions  50  pm  x  p50  m  -  2.5  10_5  cm2  the 
average  number  of  carriers  in  the  process  of  tunneling  at  any  time  is 
therefore  0.05  -  0.075,  and  this  is  indeed  much  less  than  unity. 
Nevertheless,  if  the  area  of  the  device  exceeds  5  10*4  cm2,  Eqs.  (190)  - 
(193)  require  an  additional  factor  e/tl^A  which  makes  the  noise  spectral 
density  proportional  to  It  and  A,  rather  than  to  the  square  of  these 
quantities. 

The  photoelectric  current  will  reflect  the  fluctuations  in  the  number 
of  photons  arriving  from  the  radiation  background.  The  quantum  efficiency 
will  not  exhibit  considerable  quantum  1/f  noise,  because  the  generated 
carriers  will  be  corrected  with  certainty.  Therefore  the  collection  of 
photoelectrically  generated  carriers  is  not  controlled  by  any  cross  section 
or  process  rate  affected  by  considerable  quantum  1/f  noise.  If  we  neglect 
the  1/f  noise  generated  in  the  series  resistance  of  the  diode,  there  should 
be  no  photoelectrically  generated  1/f  noise  from  a  short  -  circuited  diode. 

Since  the  various  dark  current  fluctuations  with  1/f  spectrum  are 
statistically  independent,  the  total  1/f  noise  is  simply  obtained  by 
summing  all  contributions. 


D.  1/f  NOISE  LIMITED  PERFORMANCE  OF  MIS  DIODES 

From  Eq.  (2.4)  we  write  the  total  dark  current  fluctuation  in  the  form 
Sid  *  Sid  +  Sip  +  8I5  +  8ltb  +  Slfc  +  8l{gc  +  8(qri<X>),  (1 94) 

and  the  spectral  density  of  current  fluctuations  will  be  neglecting  8(qTiO), 
S|d  ■  S,d  +  S|r  +  S|s  +  S|tb  +  Site  +  S|tsc-  (195) 
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Here  we  have  lumped  the  recombination  current  on  the  back  surface  lb 
together  with  the  surface  recombination  (generation)  current  ls.  If  we 
denote  all  the  corresponding  spectral  densities  of  fractional  fluctuations 
by  a  prime,  S'|j  -  S|j/lj2,  we  obtain 

S'ld  *  Odif^dl^S'dif  +  (*dep^d)^S'idep  +  Os^d)^S'ls 
+  (ltb/'d)2S'ltb  +  (|tc/,d)2S'ltc  +  (Wd)2S,|tsc-  (196) 

This  equation  was  obt^'rf  1  by  dividing  the  previous  equation  through  Iq2, 
and  shows  that  the  bigg^jt  contribution  will  not  necessarily  come  from 
the  process  with  the  highest  fractional  quantum  1/f  noise,  i.  e.,  with  the 
highest  1/f  noise  coefficient.  The  weight  of  each  type  of  noise  is 
determined  by  the  corresponding  squared  current  ratio. 

The  detectivity  of  infrared  detectors  is  limited  in  general  by  three 
types  of  noise:  (i)  current  noise  in  the  detector,  (ii)  noise  due  to 
background  photons  (photon  noise),  (iii)  noise  in  the  electronic  system 
following  the  detector.  We  shall  neglect  here  the  background  photon  noise 
and  the  noise  in  the  electronic  system.  The  detectivity  is  defined  as 

D*(M)  -  (AAf)1/2/NEP(cmHz1/2/W)  (197) 

where  A  is  the  area  of  the  detector,  NEP  the  noise  equivalent  power 
defined  as  the  r.m.s.  optical  signal  of  wavelength  required  to  produce 
r.m.s.  noise  voltage  (current)  equal  to  the  r.m.s.  noise  voltage  (current)  in 
a  bandwidth  Af,  and  f  is  the  frequency  of  modulation.  The  noise  equivalent 
power  NEP  is  given  by 

NEP  -  (hv/iiq)(S|d(f)Af)1/2.  (198) 

Therefore  we  obtain  for  the  detectivity 

D*(X,f)  -  (qqX/hc)[A/S|d(f)]1/2  -  (qWhc)[S|d(f)]-1/2  (199) 

We  notice  that  D*(X,f)  is  proportional  to  X  up  to  the  peak  wavelength  Xc. 

For  X  >  Xc  we  have  ti  ■  0  and  thus  D*(X,f)  -  0.  By  substituting  our  result  for 
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Sid,  we  obtain  the  general  expression  of  the  detectivity  as  a  function  of 
various  parameters  of  the  MIS  device. 

Let  us  now  evaluate  the  spectral  density  S'(f)  of  fractional 
fluctuations  in  the  various  dark  current  noise  contributions  per  cm2  of 
transversal  detector  (or  gate)  area,  including  also  a  numerical  example  for 
a  MIS  infrared  detector.  For  a  given  detector,  this  needs  to  be  divided  by 
the  area  A  of  the  detector  to  yield  the  corresponding  fractional  spectral 
densities:  S(f)  *[S'(f)]/A.  Fractional  fluctuations  are  dimensionless,  so 
S’(f)  will  have  the  dimension  of  a  reciprocal  frequency,  times  a  squared 
length  unit  which  simplifies  when  we  divide  by  the  area  of  the  detector  at 
hand.  Let  S'|d  be  the  spectral  density  of  fractional  fluctuations  in  the 
noise  caused  by  quantum  1/f  fluctuations  in  diffusion,  S'|r  in  bulk 
recombination,  S'|s  in  surface  recombination,  and  S'|t  in  tunneling.  With 
mp*  -  0.55  m0,  mn*  =  0.02  m0,  xn  -  10‘6s,  Eg  -  0.1  eV,  3kT/2  =  0.01  eV,  we 
obtain  for  a  p-type  MIS  device  with  wp»L<j 

S'ld  =*  (aHnd  +  aHnr)[e/ftnld]F(a)  -  acoh{e1/2/[f(kTpxn)1/2Np]}F(a)/a 

=(4.6  10-3/4fNp)  4  lO-10Cl/2/[10-6s  1.5  105(cm2/Vs)  4  10'21  J]1'2 
=1.8  10-6cm2/f,  [or  10'8  cm2/f,  with  ap  =2  10’5  for  incoher.  noise]  (200) 

S'ir  =  [aHee/f(*no  +  tpo)lr]tanhx  =  [otHee/(feAwn  tanhx)]tanhx 

=  otHe/fAwnj  =  4.6  10-9cm2/f  (x  »  eV/2kT)  (201) 

S’|S  =  (4a/3jc)(2/m*c2)[3kT/2  +  eU/2  +  0.1  Ve]  [e  tanhx  /f(xno+xpo)ls] 

=(4a/37c0.02)(2/500, 000) [0.025+0. 5+0. 5]  [e  tanhx/feAwnj(e*-1 )] 

=  7  10*8cm2/f  =  S|b  (202) 

S'|,b  -  (4a/3jt)(Eg  +  3kT/2)/m*c2  »  (4/9.5  137  0.02)(0.1 1/500,000) 

=  3.3  10-8cm2/f  (203) 

S’|tc  =  (4a/3jc)(Eg  +  3kT)/2m*c2  -  (4/9.5  137  0.02)(0. 1 2/1 0®) 

=  1.8  10-8cm2/f  ,  S'hsc 


(204) 
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S'|d  was  calculated  in  the  small  bias  limit  for  wp»L,  but  wp*  0.25  L  gives 
the  same  result;  the  incoherent  case  with  a  lattice  constant  of  0.65  nm 
and  0  -  320  K  was  also  listed  above  (because  a  10  |im  thick  device  is  very 
small,  so  it  may  be  applicable),  and  would  give  1.8  10*1°cm2/f  for  a  n  - 
type  device.  Eqs.  (201)  -  (204)  would  be  reduced  mp*  /mn*  *  27.5  times 
for  n  -  type  devices.  We  mention  that  S'|s  has  been  calculated  with  the 

inclusion  of  a  term  of  10%  of  the  applied  gate  voltage  V  into  the  kinetic 
energy  of  the  carriers  at  the  surface,  and  that  for  the  back  surface 
recombination  current  this  term  has  to  be  dropped  in  the  similar 
expression  of  S'|b.  However,  we  have  neglected  this  here,  because  the 

surface  recombination  terms  will  not  turn  out  to  be  important,  as  we  will 
see  below.  The  applied  gate  voltage  was  taken  to  be  V  *  5  V.  Using  Eq. 

(4.3)  and  the  current  densities  evaluated  in  Eqs.  (4.10)  -  (4.14)  to 
calculate  the  fraction  of  each  current,  we  obtain 

1  cm-2  f  S'|(f)  *  (20/132)2  1.8  10'6  +  (10/132)2  4.6  10'9  +  (3.6/132)2 
7  10-8  +(0.01/132)2  3.3  10'8  +(80/132)2  -|.8  lO’8  +(17.5/132)2  1.8  10’8 
-  3.67  10-8  +  2.6  10-H  +  5.2  10*H  +  1.9  10*16  +  6.61  10’9  +  3.17  10-10 
=  4.37  10‘8,  or  for  incoherent  1/f  noise,  7.1  10-9(p)  and  3  10*1  °(n)  (205) 

This  vaiue  can  be  used  in  order  to  estimate  the  detectivity  of  the 
device  in  our  example.  Substituting  into  Eq.  (199),  we  obtain  with  a 
quantum  efficiency  ri  =  0.7  and  wavelength  of  X  =10  pm 

D*(X,f)  =  (riqX./hc)[Sid(f)]_1/2=*  [0.7  1.6  10’19c  10*5m/(6.6  10-34js  3 

1 08m/s)]  [f/(4.37  10-8cm2  1.74  10-6A2/cm4)]1/2  ,2  107  (cm  Hz1/2/w) 
xf1/2(  [or  for  incoherent  1/f  noise,  5  107  (p),  and  2.5  108  (n)].  (206) 

In  conclusion  we  note  that  for  the  relatively  large  devices  which  we 
have  considered,  most  of  the  quantum  1/f  noise  co.nes  from  fluctuations 
in  diffusion  and  in  the  tunneling  rate  via  impurity  centers  in  the  band  gap. 
The  effective  mass  of  the  carriers  is  present  in  the  denominator  of  all 
quantum  1/f  noise  contributions  except  the  coherent  quantum  1/f 
fluctuation  present  in  the  diffusion  current  of  large  devices.  In  smaller 
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devices  the  diffusion  current  will  also  be  given  by  the  conventional 
quantum  1/f  formula  which  contains  the  effective  mass  of  the  carriers  in 
the  denominator.  For  Umklapp  scattering  the  mass  of  the  carriers  in  the 
denominator  is  even  squared.  Consequently  we  expect  lower  quantum  1/f 
noise  from  n  -  type  devices,  in  which  the  minority  carriers  are  holes, 
particularly  if  the  devices  are  very  small,  e.g.,  below  10pm. 

E.  DISCUSSION  AND  RECOMMENDATIONS 

The  transition  from  coherent  state  quantum  1/f  noise  to  conventional 
quantum  1/f  noise  is  particularly  interesting,  and  should  be  studied 
experimentally.  This  is  possible  with  a  sequence  of  devices  of  smaller 
and  smaller  size,  and  will  show  a  considerable  change  in  noise  at  a  size  of 
the  order  of  10  m2.  The  theory  of  the  transition  is  not  yet  well  developed. 
Therefore,  this  experiment  nas  particular  importance;  we  do  not  know  if 
the  parameter  s  is  sufficient  to  characterize  the  transition,  and  if  the 
parameter  s  should  not  be  replaced  by  a  power  of  s,  or  by  any  other 
function  of  s.  The  interpolation  formula  used  here  is  just  a  guess,  or  a 
speculation  guided  by  the  physical  understanding  of  coherent  quantum  1/f 
noise  as  a  collective-field  effect,  and  of  conventional  quantum  1/f  noise 
as  an  effect  which  is  not  based  on  the  collective  field  state  of  the 
particles,  but  arises  from  the  individual  field  of  each  carrier. 

The  most  interesting  component  of  the  recombination  current  is  the 
surface  recombination  current  which  plays  a  major  role  in  the  case  of 
infrared  detectors  with  pn  junctions.  In  the  case  of  MIS  devices  this  role 
is  not  so  important,  as  our  calculation  shows.  Nevertheless,  one  should 
try  to  reduce  both  the  concentration  of  recombination  centers  and  the 
value  of  the  surface  potential  jump  U.  This  can  be  accomplished  with 
careful  surface  treatment,  and  with  a  good  passivation  layer.  Si02  layers 

have  been  successfully  used  by  Radford  and  Jones  in  ion-implanted  and 
double  -  layer  epitaxial  HgCdTe  photodiodes34. 

In  general  the  larger  life  time  of  the  carriers  in  MIS  devices,  compared 
to  junction  devices  is  due  to  the  absence  of  the  damage  inflicted  by  ion  - 
implantation,  or  by  the  heavy  doping  required  in  double  -  layer  epitaxial 
photodiodes.  The  quantum  1/f  noise  is  inversely  proportional  to  this  life 
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time.  Therefore,  MIS  devices  should  have  lower  1/f  noise.  On  the  other 
hand,  1/f  noise  present  in  the  applied  gate  voltage,  in  the  timing  of  the 
readout  and  the  value  of  the  readout  potential  will  be  added  as  a  1/f  noise 
source,  if  it  is  present.  In  the  present  calculation,  however,  this  noise 
source  has  not  been  included. 

Any  reduction  in  the  concentration  of  tunneling  centers  present  in  the 
band  gap  will  have  a  positive  effect  on  quantum  1/f  noise.  As  we  have 
seen  in  Sec.  III.4B,  p  -  type  devices  should  yield  less  tunneling  via  band 
gap  centers.  The  effective  mass  present  in  the  denominator  of  the 
quantum  1/f  noise  formula  in  this  case  should  just  be  the  effective  mass 
of  the  carriers  after  the  tunneling  process,  i.e.,  the  effective  mass  of  the 
outgoing  carriers  emerging  from  the  process  we  have  considered,  or  the 
effective  mass  of  the  carriers  coming  in  to  the  process  of  tunneling 
toward  the  centers  in  the  band  gap.  Here  we  have  considered  the  tunneling 
process  as  the  slower  process  which  actually  controls  the  rate  of 
tunneling  via  band  gap  states.  The  capture  of  carriers  by  the  band  gap 
states  is  the  second  part  of  this  compound  process  and  has  been 
considered  fast  enough,  so  that  it  does  not  limit  the  rate  of  the  total 
process.  In  general,  however,  both  parts  of  the  process  have  to  be 
considered  as  a  limitation  on  the  rate,  and  in  this  case  our  noise  formulae 
have  to  be  revised  through  the  inclusion  of  an  additional  term  similar  to 
the  recombination  noise  term. 

In  the  case  of  very  small  MIS  devices,  where  only  conventional  quantum 
1/f  noise  should  be  present,  we  may  find  lower  noise  in  the  n  -  type 
devices,  whose  bulk  minority  carriers  are  holes  with  much  larger 
effective  masses  than  the  electrons.  This  may  happen  in  spite  of  the 
larger  tunneling  via  band  gap  centers  located  right  under  the  surface  of 
these  devices. 

Finally,  we  would  like  to  emphasize  that  the  present  study  has 
attempted  to  explain  the  basic  concepts  of  quantum  1/f  noise  and  to 
illustrate  their  application  to  MIS  infrared  detectors.  Although  we  have 
tried  to  pursue  the  calculation  all  the  way  to  the  evaluation  of  the 
detectivity,  the  data  which  we  used  in  the  calculation  may  not  be 
applicable  in  the  practical  case  at  hand,  and  may  have  to  be  replaced  with 
pertinent  data  in  any  concrete  case. 
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111.5.  QUANTUM  1/f  NOISE  IN  SQUIDS 


As  we  have  seen  in  Sec.  11.1-3  above,  any  cross  section  or  process 
rate  defined  for  electrically  charged  particles  must  fluctuate  in  time 
with  a  1/f  spectral  density  according  to  quantum  electrodynamics,  as  a 
consequence  of  infrared-divergent  coupling  to  low-frequency  photons. 

This  fundamental  effect  leads  to  quantum  1/f  noise  observed  in  many 
systems  with  a  small  number  of  carriers,  and  is  also  present  in  the  cross 
sections  and  process  rates  which  determine  the  resistance  and  tunneling 
rate  in  Josephson  junctions,  providing  a  lower  limit  of  the  observed  1/f 
noise. 

In  a  Josephson  junction  the  normal  resistance  Rn  of  the  barrier  is 
proportional  to  a  scattering  cross  section  or  transition  rate  experienced 
by  the  electron  in  quasiparticle  tunneling  and  by  the  Cooper  pairs  below 
the  critical  current  !c.  Therefore 

Rrr2SRn(f)  =  (4<*/3t0  [(Av)2/C2Nf) 

=  (8u/3k)  (vF2/c2Nf)  =  4  10-14/fn  (207) 

where  we  have  approximated  (Av)2  with  2vF2,  vF  being  the  Fermi  velocity, 
and  the  number  of  carriers  N  simultaneously  present  in  the  barrier  of 
volume  Q  (in  p3)  by  1 07,  for  barriers  wider  than  10-7cm. 

Assuming  a  linear  relationship  between  the  critical  current  lc  and 
Gn^n'1.  we  obtain  similar  to  Rogers  and  Buhrman4^,  substituting 
however  our  quantum  1/f  source  given  by  Eq.  (207)  for  Rn,  the  spectral 
density  of  voltage  fluctuations 


S v(f)  -  (4/f)  10-12(T/3K)[Rsg(V)/(Rs+Rj)]2 
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*[lcRn(l2/lc2  -1)-1/2+  g(V)V]2Q-l,  (208) 

where  Rj(V)  is  the  junction  resistance,  Rs  the  shunt  resistance,  and 
g(V)*Rn/Rj. 

The  noise  caused  in  a  SQUID  by  the  source  considered  above  can  be 
obtained  as  the  sum  of  the  noise  contributions  from  the  two  junctions. 

The  above  quantum  1/f  results  of  Eqs.(207)  and  (208)  are  in  good 
quantitative  agreement  with  the  experimental  data. 

In  conclusion,  the  fundamental  quantum  1/f  fluctuations  of  the 
cross  sections  and  transition  rates  which  determine  the  normal 
resistance  have  been  evaluated  in  this  Section  for  the  case  of  a  Josephson 
junction.  Considering  the  velocity  change  in  the  quantum  1/f  formula 
equal  to  twice  the  Fermi  velocity  and  the  concentration  of  carriers  in  the 
barrier  1019cm'3,  a  spectral  density  of  fractional  fluctuations  in  the 
normal  resistance  of  the  barrier  of  4  10'14/f  was  obtained  for  a 
Josephson  junction  with  a  volume  of  the  barrier  of  10-12cm3.  These 
fluctuations  are  inversely  proportional  to  the  barrier  volume  and  result  in 
voltage  fluctuations  both  directly  and  through  the  dependence  of  the 
critical  current  on  the  normal  resistance,  in  good  agreement  with  the 
experimental  data. 

IV.  SUMMARY  OF  OPTIMAL  DESIGN  PRINCIPLES 


As  a  general  conclusion  of  the  present  study,  the  following  general 
principles  of  optimal  quantum  1/f  noise  reduction  emerge: 

1.  Avoidance  of  coherent  state  quantum  1/f  noise  by  device  size 
reduction  below  the  coherent  limit.  This  size  limit  is  concentration- 
dependent,  as  seen  from  the  expression  of  the  coherence  parameter  s  * 
2e2N'/mc2  -  5  10'13  cm^x  N'  defined  in  Eq.  (80).  N’  *  nA  is  the  number  of 
carriers  per  unit  length  of  the  device  in  the  direction  of  current  flow.  A 
is  the  cross-sectional  area  of  the  current-carrying  device,  and  n  is  the 
concentration  of  carriers.  For  s«1  we  expect  conventional  quantum  1/f 
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noise,  while  for  s»1  the  much  larger  coherent  state  quantum  1/f  noise  is 
to  be  expected. 

2.  Avoid  control  of  the  device  current  or  voltage  by  elementary  cross 
sections  or  process  rates  tested  by  a  small  number  of  carriers  only. 

Indeed,  the  number  of  carriers  interrogating  the  cross  section  or  process 
rate  appears  in  the  denominator  of  both  the  conventional  and  coherent 
quantum  1/f  noise  formulae.  In  particular  avoid  current  concentrations  in 
bottlenecks,  and  current  inhomogeneities.  In  junction  devices  higher 
lifetimes  of  the  carriers  lead  to  an  increase  in  the  number  of  carriers 
present  in  the  sample  which  have  tested  the  current-controlling  cross 
sections,  and  therefore  lead  to  lower  quantum  1/f  noise. 

3.  Avoid  control  of  a  device  exhibiting  conventional  quantum  1/f 
noise  through  elementary  processes  which  involve  large  accelerations  of 
the  current  carriers,  or  large  velocity  changes.  The  squared  vector 
velocity  change  appears  as  a  factor  in  the  conventional  quantum  1/f  noise 
formula.  For  example,  Umklapp  scattering,  inter-valley  and  lattice 
scattering  are  respectively  worst,  very  bad,  and  bad,  compared  with 
ionized  impurity  scattering,  in  terms  of  the  fractional  mobility 
fluctuations  which  they  yield.  For  a  given  scattering  mechanism,  choosing 
current  carriers  with  a  large  effective  mass  will  in  general  reduce  the 
conventional  quantum  1/f  noise,  because  for  the  same  momentum 
transfers  this  results  in  smaller  accelerations.  Bulk  recombination 
control  of  the  current  through  a  pn  junction  will  lead  to  lower  quantum 

1/f  noise  than  having  the  current  controlled  even  in  part  by  surface 
recombination,  because  the  surface  recombination  centers  are  in  a  high 
localized  field  3gion  at  the  interface  between  bulk  and  the  passivation 
layer.  Therefore,  the  best  passivation  is  one  which  reduces  the  number 
and  the  cross  section  of  the  surface  recombination  centers,  while  also 
providing  the  smallest  surface  potential  jump. 

Consequent  use  of  these  principles,  leads  to  lower  1/f  device  noise. 
The  quantum  1/f  theory  allows  for  CAD  optimization  of  1/f  device  noise 
suppression. 
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